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We discuss physical properties of 'integer' topological phases of bosons in D=3+l dimensions, 
protected by internal symmetries like time reversal and/or charge conservation. These phases invoke 
interactions in a fundamental way but do not possess topological order and are bosonic analogs of free 
fermion topological insulators and superconductors. While a formal cohomology based classification 
of such states was recently discovered, their physical properties remain mysterious. Here we develop 
a field theoretic description of several of these states and show that they possess unusual surface 
states, which if gapped, must either break the underlying symmetry, or develop topological order. In 
the latter case, symmetries are implemented in a way that is forbidden in a strictly two dimensional 
theory. While this is the usual fate of the surface states, exotic gapless states can also be realized. 
For example, tuning parameters can naturally lead to a deconfined quantum critical point or, in 
other situations, a fully symmetric vortex metal phase. We discuss cases where the topological 
phases are characterized by quantized magnetoelectric response 8, which, somewhat surprisingly, 
is an odd multiple of 2n. Two different surface theories are shown to capture these phenomena - 
the first is a nonlinear sigma model with a topological term. The second invokes vortices on the 
surface that transform under a projective representation of the symmetry group. A bulk field theory 
consistent with these properties is identified, which is a multicomponent BF theory supplemented, 
crucially, with a topological term. A possible topological phase characterized by the thermal analog 
of the magnetoelectric effect is also discussed. 
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I. INTRODUCTION 

In the last several years we have learnt of many phases 
and phase transitions of quantum many body systems 
that go beyond conventional concepts of broken symme- 
try Striking examples are gapped phases of matter with 
topological order. These are characterized by emergent 
excitations with unusual quantum statistics and ground 
state degeneracies that depend on the topology of the un- 
derlying manifold 1 . Other examples are gapless phases 
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of matter where the gaplessness is protected but not by 
a broken symmetry. The most familiar example of such 
a phase is a Fermi liquid but gapless spin liquids and 
various non-fermi liquid phases provide other examples. 
Other examples of exotic gapless states include decon- 
fined quantum critical points^, which nevertheless may 
occur between conventional ordered phases. It has be- 
come clear that all these "exotic" phases have some emer- 
gent non-local structure in their ground state wavefunc- 
tion. This non-locality is loosely captured by saying that 
the local degrees of freedom have long range entangle- 
ment in the ground state wavefunction. 

More recently it has become clear that some of the 
interesting properties of such exotic phases are also 
shared by phases with only short ranged entanglement. 
The best example is provided by free fermion insula- 
tors/superconductors with topological band structure^. 
As is well known these insulators have non-trivial pro- 
tected surface states but only have short range entangle- 
ment. 

Let us now generalize the concept of a topological 
insulator. Consider a system of interacting bosons or 
fermions in d dimensions, possibly including global sym- 
metries. Assume that there is a unique ground state 
with gapped excitations on all closed manifolds i.e. there 
is no topological order . We refer to such phases as 
Short Range Entangled ( SRE) statesP. Do distinct SRE 
phases exist that share the same symmetry, but differ at 
the level of topology? A possible distinction, for example, 
is presence of protected states at the boundaries. 

It is interesting to pose this question for interaction 
dominated phases of matter, i.e phases that require in- 
teractions between the underlying particles to stabilize 
them. A classic example (one which is the subject of this 
paper) are gapped phases of interacting bosons. When a 
global U(l) symmetry is present, associated with boson 
number conservation, such phases correspond to bosonic 
insulators. For bosons the non-interacting limit is a sim- 
ple condensate so that interactions are necessary to sta- 
bilize insulating phases. 

Recently there has been much progress in identify- 
ing such topological phases of bosons in diverse dimen- 
sions. In one spatial dimension (ID), these include 
states like the Haldane (or AKLT) state of gapped spin- 
1 chains^. Using the matrix product representation of 
gapped states^^ they are argued to be classified by pro- 
jective representations of the symmetry group (G) or 
equivalently by the second group cohomology of sym- 
metry group G. While in 2D and 3D such rigorous re- 
sults are not available, recently, Chen, Liu, Gu and Wen 9 
have proposed that SRE topological phases of bosons are 
captured by the higher dimensional cohomology groups. 
While Chen et alP restrict attention to the non-chiral 
subset of these states (i.e. ones that do not have a net 
imbalance of left and right movers at the edge of ) pro- 
tected by symmetry, Kitae\ff^ has also considered chiral 
states. An explicit example of such a phase in d = 2 was 
constructed by Levin and GrPM Later, in two dimensions 




FIG. 1: Schematic depiction of a 3D symmetry protected 
topological phase, with two conserved species of bosons 
(U(l)xU(l)) and time reversal symmetry (Zj). The bulk 
is insulating and corresponds to a condensate of vortices of 
both species (shown as black and white loops). In the topo- 
logical phase, the vortex line of one species that ends on a 
surface carries half charge of the other species. Such a sur- 
face, it may be argued, does not have a trivial gapped phase, 
where the symmetries are preserved. Pictorially, the vortex 
lines may be viewed in ID Haldane phases, with half charged 
end states. 



many of these results were obtained in a simpler way us- 
ing a Chern-Simons approaclf^. The classification shows 
that interaction-dominated phases of bosons with a bulk 
gap and a unique ground state on closed manifolds exist 
that have non-trivial surface states. For instance in d = 2 
bosons with a conserved global U(l) but no other sym- 
metries have an "integer quantum Hall state" with non- 
trivial protected edge states but no fractional bulk excita- 
tions. In particular, the quantized Hall conductance was 
predicted^ to be an even integer in units of q 2 /h, where 
q is the elementary charge of the bosons. Very recently 
a simple physical realization of such a bosonic integer 
quantum Hall effect has been provided, and its proper- 
ties shown to agree with the formal classification^. In 
general, depending on the detailed symmetry group and 
dimensionality of the system, a variety of such "Sym- 
metry Protected Topological" (SPT) phases have been 
argued to exist. In d = 1 or in d = 2 we have a fair un- 
derstanding of the universal physical properties of these 
SPT phases. However in d = 3, though the cohomology 
classification allows for the existence of a number of dif- 
ferent SPT phases, including bosonic generalizations of 
topological insulators, their physical properties have not 
thus far been elucidated. This will be our primary task. 
Potentially, the approach here could also serve to inde- 
pendently classify these phases, which however is left for 
future work. 

Surface States and Projective Vortices: Typically, 
physical manifestations of topological phases involve pro- 
tected surface states and/or a quantized response. We 
show that 3D SRE bosonic topological phases can also 
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be characterized in this way. First, we will see that such 
phases only occur in the presence of symmetry. Then, 
as with free fermion topological insulators, the surface 
states are either gapless or, if gapped, must develop ei- 
ther conventional or topological order. Conventional or- 
der involves breaking the underlying symmetry defining 
the phase. The other possible gapped state involves topo- 
logical order that develops just at the boundary. While 
we discuss multiple routes that reach this conclusion, a 
particularly intuitive picture is the following. Consider 
for a moment a system with two independent species 
of bosons (species 1,2 whose numbers are independently 
conserved by U(l)xU(l) symmetry) as well as time re- 
versal (Z2) symmetry. We describe a three dimensional 
topological insulator of such a system where the surface 
has the following unusual property. If bosons of species 1 
is condensed on the surface, then its vortices carry half- 
charge of species 2, and vice versa. This implies that 
disordering the surface superfluid, which corresponds to 
condensing vortices, also condenses species 2. That is, 
there is an obstruction to realizing a fully symmetric and 
gapped surface. It is important that the vortices carry 
fractional charge of species 2, which cannot be screened 
by bosons on the surface. Time reversal plays an impor- 
tant role in fixing this charge to be ±1/2. Some of the 
options then are (i) to break U(l) symmetry on the sur- 
face such as in the surface superfluid or (ii) to condense 
pairs of vortices that carry no charge, which leads to 
Z2 topological ordei^31, or (iii) both vortices and bosons 
are gapless but not condensed. For bosons, which are 
typically either gapped or condensed, this corresponds 
to a (multi) critical point. In fact, it corresponds to a 
deconfined critical point, analogous to that proposed in 
the context of a continuous quantum phase transition be- 
tween a Neel and valence bond solid (VBS) on the square 
lattice^. There it was pointed out that vortices of the 
VBS order carry spin 1/2, i.e. 'half charge' of the other 
symmetry^. More exotic gapless Bose or vortex metal 
phases are also possible. The key point however is that a 
trivial gapped surface state preserving all symmetries is 
not possible. A gapped surface that preserves all sym- 
metries, necessarily exhibits topological order. Moreover, 
this topologically ordered state differs from a truly two 
dimensional phase with the same topological order and 
symmetry. We discuss specific examples where the ac- 
tion of symmetry on the surface topological excitations 
is distinct from that of any two dimensional state. 

Finally, a different symmetry breaking option is to 
break time reversal, which can lead to a gapped surface. 
Given that this is an insulator, with gapped bosons both 
in the bulk and on the surface, it can be characterized by 
a quantized response - the magneto electric polarizability 
9. This survives on reducing to a single conserved charge 
(a charge U(l) along with time reversal symmetry) where 
this leads to 9 = 2n, which is a defining property of this 
phase as elaborated below. 

Magnetoelectric polarizability of bosonic insulators: It 
is well known that free fermion topological insulators, 



with broken time reversal on the surface that removes the 
surface states, have a quantized magneto-electric effeclP^ 
captured by the topological theta term: 

Co = ^E.B (1) 

(we have set h = c = e = 1 and E, B are applied electric 
and magnetic fields). For fermion topological insulators 
9 = 7r (mod27r), corresponding to a half integer Hall effect 
on the surface. The 2-7T ambiguity in 9 corresponds to the 
fact that one may deposit an integer Quantum Hall layer 
on the surface^. 

One of our central observations is that for 3D bosonic 
topological insulators, with conserved charge and time re- 
versal symmetry in the bulk, 9 is only defined modulo 
4-7T, and the topological phase corresponds to 9 = 2tt. 
This implies, for example, the domain wall between op- 
posite time reversal symmetry breaking regions on the 
surface induces a protected mode, the edge state corre- 
sponding to the o~ X y ~ 2q 2 /h quantized Hall effect of 
bosons. The 4-7T ambiguity in 9 corresponds to the fact 
that one may deposit an integer Quantum Hall layer of 
bosons on the surface, which must have an even integer 
Hall conductance. The gapless surface state may there- 
fore be viewed as a Quantum Hall state that fluctuates 
between o~ xy = ±q 2 /h. A powerful route to obtaining a 
theory of such a state is to start with the theory of the 
quantum phase transition between distinct integer quan- 
tum Hall states of bosons. Such an approach when ap- 
plied to fermionic topological insulators correctly yields 
the single Dirac cone surface state. We therefore con- 
struct a network model to describe the bosonic integer 
quantum Hall transition in two dimensions and use that 
to derive a field theory of this transition. This then nat- 
urally leads us to the right effective field theory of the 
surface of the three dimensional bosonic topological in- 
sulator. We show that that key properties, such as the 
presence of projective vortices, emerge from this starting 
point. 

Possible Phases with Half Quantized Surface Thermal 
Hall effect: Note, in the discussion above, the quantized 
Hall effect of d = 2 SRE bosons immediately constrained 
the physics in 3D. An analogous argument for thermal 
Hall conductance can be made, which in d = 2 is quan- 
tized to 8 times the quantum of thermal conductance: 

7T 2 k 2 

Kxy/T — 8n 3h B . A realization of n = 1 is the Kitaev Eg 
state, with 8 chiral bosons at the edge. Therefore, a three 
dimensional phase protected by time reversal symmetry 
can be conceived, that has surface states at a domain wall 
between opposite T breaking regions on the surface, with 
the 8 chiral boson modes of the E% state. Indeed we will 
construct field theories of such states - whether these can 
be realized in microscopic models regularized on a lattice 
is an open question that we will briefly discuss. We note 
that if indeed this particular class of states exists then it 
is likely to be outside the cohomology classification. 

3D 'BF+FF' Field Theory: The general arguments 
above will be shown to be consistent with the following 
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d = 3 field theory: 

2nC 3D = ^ Xa B^d x ai + 6 £ ^e^^/^J 
i i, J 

where the index I refers to boson species, and bosons 
four currents are represented by j M/ = ^ze^ uXo d v B\^ a 
while the curl of a 1 represents the vortex lines. The first 
'BF' terrrPD j us t represents the 2ir phase factor of tak- 
ing a particle around its vortex. The key topological 
properties however are determined by the second term, 
which attaches quantum numbers to vortices . To avoid 
topological order at the surface, and ensure bosonic ex- 
citations, detK — 1 and diagonal entries of K are even 
integers. In most cases we will take two species with 
K = a x . Furthermore, time reversal symmetry con- 
strains = 0, 7r, the latter being the topological phase. 
This 9 for the internal gauge fields a 1 associated with 
vortex lines should be distinguished from the 6 for the 
external electromagnetic field we discussed above. Cou- 
pling to an external electromagnetic field allows one to 
obtain the quantized magneto electric effect discussed 
above. Related theories have appeared in the context 
of three dimensional topologically ordered phases and 
superconductors^ where only the first BF term in EqnJS] 
appears with a different coefficient. 

Recently, it was proposed that free fermion topological 
insulators are captured by similar theories,^ again with 
only the first term of Eqn. [2] and coupling to the external 
field. The surface states in su ch th eories were argued 
to be bosonized Dirac fermiond 20 * 21 *. While this is an 
intriguing idea, we point out certain problems with the 
specific implementation in Ref. [20l In particular, we find 
no evidence for surface fermions (see Appendix [A]) in the 
single component BF theory of Ref. 1201211 It is presently 
unclear if there is a suitable modification that overcomes 
these issues. Finding an effective field theory description 
for 3D fermionic topological insulators remains an open 
problem. On the other hand the field theory discussed 
in Ref. [22] retains only the second term, which leads to 
gapless excitations in the bulk, as noted in Ref. [33J which 
differ from the gapped phases of interest here. Thus, it 
will be important to combine both these terms. Further 
we will need to allow for multiple species of bosons. This 
is reminiscent of the K-matrix description of several SPT 
phases in 2Eff2l, where a two component Chern-Simons 
theory, arising from two species of bosons, was required 
to capture topological properties. 

We also display continuum field theoretic models in 
D = 3 + 1 dimension that realize some of the topological 
phases we describe. These are obtained as perturbations 
of non-linear sigma models in the presence of a topologi- 
cal theta term. The theta term has the effect of endowing 
topological defect configurations of the continuum fields 
with non-trivial global quantum numbers. We show the 
connection to the topological BF theory and to the the- 
ory of protected surface states. 



II. OVERVIEW 

As this paper is long and discusses bosonic SPT phases 
from several points of view it is helpful to provide an 
overview. Let us first recall the defining properties of 
such a phase. The first requirement is that it has a 
bulk gap but no "intrinsic topological order" implying 
a unique ground state on all closed manifolds (loosely 
characterized as SRE). The second requirement is that 
in the presence of a boundary the surface cannot be in 
a trivial gapped insulating phase that preserves all the 
symmetries. One of our goals in this paper is to seek 
an effective "Landau-Ginzburg" description of the sur- 
face theory which naturally incorporates this restriction. 
Below we motivate our approaches to obtaining such an 
effective theory. 

It is worthwhile to emphasize a different point of view 
on SPT phases. One definition of SRE is that the ground 
state wavefunction can be smoothly deformed into a 
product of local degrees of freedom without undergoing a 
phase transition. (This definition - used by Chen et alP- 
is a bit more restrictive than an alternate one used by Ki- 
taev). An SPT phase may then be described as a phase 
with a ground state wavefunction which can be smoothly 
deformed into a product state if and only if the defining 
global symmetry is broken. Let us now briefly discuss 
the surface states in the presence of a boundary from 
this point of view. If the surface cannot be in a trivial 
gapped state then what else can it be? Clearly the sur- 
face may spontaneously break the symmetry. Then in the 
absence of the symmetry in the bulk there is nothing spe- 
cial about the surface consistent with it being deformable 
into a product state. For a d — 3 system the surface is 
two dimensional and there is the interesting possibility 
that the surface itself has intrinsic topological (or other 
more exotic) order, for instance, a fractionalized insula- 
tor described by a 2 + 1 dimensional deconfined Zi gauge 
theory. Can we envisage SPT phases where the surface 
necessarily has such intrinsic topological order though the 
bulk does not? The answer is no. If we add a symmetry 
breaking perturbation to the bulk we should be able to 
get a trivial surface for an SPT phase. So the surface 
theory should be such that it allows for intrinsic surface 
topological order but does not require it. The effective 
field theories we obtain for the surface state will all have 
these general features. 

To motivate our line of attack on 3D SPT phases we 
adopt in this paper, it is helpful to recall the physics in 
lower dimensions, particularly in 2D. There, a number 
of different symmetry protected topological phas es can 
be represented by the same 2 x 2 K matrix^!!] K = 
o x , and differ only in the way the edge fields transform 
under different symmetries. The K matrix imposes the 
following term in the edge Lagrangian: 

£cd g ciD = ^-d t 4>id x (f)2 - V u d x (j)id x (t)j 

Z7T 

the first term is determined by the bulk K matrix and 
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fixes the commutation relations between the edge fields 
01,2- The second term is non- universal and depends on 
the detailed interactions. This is the Lagrangian for a 
1+1D Luttinger liquid of bosons, where the two phases 
4>i may, for instance, be identified with the boson phase 
(<p) and twice of the conjugate phase {29) related to den- 
sity wave order. However, in the edge theory above, the 
transformation property of the fields is determined by 
the topology of the bulk. Thus, f or exa mple in the inte- 
ger quantum Hall effect of bosonJ^EEH ) both fields <f>\ t 2 
are translated by the U(l) charge symmetry, which does 
not occur in a purely ID Luttinger liquid where only 
one of the fields (the boson phase) is charged. Depend- 
ing on the transformation laws, cosine terms compatible 
with the symmetries may be added to the theory, which 
determines its ultimate state. In the previous example, 
with U(l) charge conservation symmetry which cannot 
be spontaneously broken in ID, the edge remains a gap- 
less Luttinger liquid. In other cases, the edge m ay be 
gapped after spontaneously breaking the symmetr ^ 11 ! 12 !. 
Nevertheless, it is useful to begin with the Luttinger liq- 
uid edge theory above as a parent theory of the edge, 
to which various anisotropics are added. Implicit in this 
approach to 2d SPT phases is that it is a useful device 
to first consider the surface theory for a system with en- 
larged U{1) x U(l) (described by the 2x2 K-matrix al- 
luded to above) , and then add anisotropy terms to break 
the symmetry to the one under consideration. 

This suggests that in three dimensions as well, we 
should first construct a model of the surface. We do 
this in two different ways. First, we argue for a parent 
theory of the surface state, analogous to the Luttinger 
liquid in one lower dimension, and determine transfor- 
mation of various fields under the symmetry operations. 
As in two dimensions we initially consider an enlarged 
continuous global symmetry U(l) x ?7(1) and eventually 
break it down to the symmetry of interest. In its minimal 
form, this leads to a deconfined quantum critical action: 



£cd g c2D = ^ K^m ~ «*/i)^v| 2 + K(dfj,a v - d v a^f 

CT=± 

where the boson density is represented as the flux of the 
gauge field n — {d x a y — d y a x )/2iT and the vortices of this 
boson field ip a are coupled minimally to the gauge field. 
The index a refers to the multiple components of the vor- 
tex field which is forced upon us due to the symmetries of 
the problem. The key observation is that since the vortex 
is not a local object, they can transform by an additional 
phase when acted by a symmetry. Such a projective rep- 
resentation of a symmetry group ensures that a trivial 
representation of the symmetry does not arise on combin- 
ing with any number of bosons. For example, projective 
representations can lead to half charged vortices that can- 
not be screened by integer charge bosons. Furthermore, 
projective representations enforce degeneracies which are 
crucial. If a single bosonic vortex species was present at 
low energies, we could perform a duality back to the usual 



boson phase variables. However, having multiple vortex 
fields demands a dual description like the one above, in 
terms of vortices coupled to a gauge field. Operators that 
insert bosons act as monopole insertion operators. In the 
context of 2D deconfined quantum criticality, the action 
above appeared where the density n ~ S z is the spin 
density of an easy plane antiferromagnet, and the gauge 
invariant combination of the vortex fields ijj+ip- ~ elS 
is the valence bond solid order parameter. However, as 
a theory of the surface on a 3D topological phase, the 
fields will transform under internal symmetries, which 
will dictate additional terms that can be added to this 
Lagrangian. Nevertheless, this will prove a useful par- 
ent theory of the 3D surface states, much like the 2D 
Luttinger liquid. 

An alternate derivation of the same surface theory is 
arrived at by recalling that the surface state of the 3D 
fermionic topological insulators, the single Dirac cone, is 
also the dispersion at the critical point between integer 
quantum Hall plateaus^ in a clean two dimensional sys- 
tem. The time reversal symmetry present in the former 
case automatically tunes the system to the criticality. By 
analogy, here we consider a network model, poised at the 
transition point between two bosonic integer Quantum 
Hall states. The resulting theory is an 0(4) non-linear 
sigma model with a topological term as in the Euclidean 
Lagrangian below: 



edge2D 



1 
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i^^ X n&d,g)(gld„g)(gld,g)} 



where the 0(4) vector has been written in terms of an 
SU(2) matrix g. Reassuringly, this has been argued in 
Ref. [2E]to be equivalent to the deconfined quantum crit- 
ical point surface theory once appropriate anisotropics 
are introduced, as we also discuss below. 

Fermionic vortices and vortex metals: A natural case 
to consider is when there is a single U(l) symmetry lead- 
ing to a conserved 'charge' and time reversal symmetry. 
Within our formalism, it is convenient to begin with two 
separately conserved species of bosons and then break 
this symmetry down to a single physical U(l) by allow- 
ing for mixing between boson species. Assuming inde- 
pendent species, the topological phase is characterized by 
vortices of one species carrying half charge of the other 
species (see figure [I]). However, if the two species are 
the same, then the resulting vortices are conveniently 
viewed as fermionic (as argued in section IV D I. Thus, 



if the surface is superfluid (while the bulk remains in- 
sulating), vortices may be thought of as point particles, 
and expected to be fermions. In the presence of a finite 
background density of charge, it may be operationally 
difficult to define fermonic statistics^!. However, if the 
conserved charge refers to S z component of spin, then 
time reversal forbids background density, which allow for 
statistics to be sharply defined. In this case the vortex 
density is time reversal invariant and generically will be 
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non-zero in the ground state. Destroying the XY order 
by proliferating the vortices will then lead to the forma- 
tion of a Fermi surface of these vortices which as usual 
will be coupled to a fluctuating U(l) gauge field. Thus 
in this case a rather exotic symmetry preserving state is 
possible: a vortex Bose 'metal^. The vortex metal is a 
compressible metallic state of bosons that preserves all 
physical symmetries. Such a state was proposed, and its 
properties analysed, in Ref. 28, as a possible quantum 
vortex liquid state in two space dimensions. Here we see 
that this rather exotic state arises very naturally at the 
surface of a 3d boson SPT phase when the surface su- 
perfluidity is destroyed. We note here that the general 
properties of doped Conformal Field Theories in 2 + 1 
dimensions and the possibility of obtaining Bose metal 
phases has also been recently discussed^. 

Surface topological order and symmetry: It will also 
be extremely useful to consider the properties of sur- 
face states which have intrinsic topological order even 
though the bulk does not. A simple but powerful exam- 
ple is a state with Z^ topological order. In cases with a 
global 17(1) symmetry these may be obtained by pairing 
and condensing vortices. A key feature of Z 2 topolog- 
ical order is the existence of three distinct topological 
quasiparticles (often denoted the electric, magnetic, and 
their composite particles). These three all have mutual 
7r statistics. Further in the common kind of Z^ topologi- 
cal order (which is what is pertinent for this paper) two 
of these are bosons while their bound state is a fermion. 
Usually a trivial phase is obtained from such a topologi- 
cal ordered phase by confinement. This occurs when one 
of the bosonic topological quasiparticles condense. When 
the topologically ordered state is obtained at the surface 
of a boson SPT phase, we will show that both types of 
bosonic topological quasiparticles transform nontrivially 
under the defining symmetry. When such a quasiparti- 
cle condenses, the surface topological order is indeed lost 
but the resulting non-topological phase breaks the global 
symmetry. 

Two dimensional bosons on clean lattices at commen- 
surate filling also exhibit Z 2 topological ordered phases. 
There again there are two bosonic topological quasiparti- 
cles - one (dubbed the chargon) which carries fractional 
boson charge and the other (usually called the vison) 
which transforms projectively under lattice symmetries. 
In the presence of these lattice symmetries, again a triv- 
ial gapped symmetry preserving insulator is forbidden as 
condensation of either chargons or visons breaks global 
symmetry. In contrast as outlined above at the surface of 
the boson SPT phase both bosonic quasiparticles trans- 
form projectively under internal symmetries and these 
are therefore sufficient to preclude the trivial gapped in- 
sulator. 

Strictly two dimensional time reversal invariant 
abelian gapped phases were recently classified in Ref. [24] 
In the specific case of Zi topological order, the realiza- 
tion of symmetry at the surface of the 3d SPT phases we 
discuss is beyond this classification. Thus the realization 



of symmetry at the surface topologically ordered state 
is distinct from what is allowed in a strictly two dimen- 
sional state with the same topological order and internal 
symmetry. These issues are discussed in Appendix [D) 

Topological Surface States and Higher Dimensional 
Lieb-Schultz-Mattis Theorems: The feature that the two 
dimensional surface state of a 3d boson topological in- 
sulator cannot be in a trivial gapped insulating phase is 
reminiscent of the Li eb-Sh ultz-Mattis (LSM) theorem^ 
and its generalization^^ to states of bosonic systems 
at a fractional filling on clean 2d lattices. Indeed in 
both cases either a symmetry must be broken or there 
is topological order, or more exotic long range entan- 
glement possibly with protected gapless excitations. It 
is thus perhaps no surprise that the field theoretic de- 
scriptions we obtain of the surface states of the 3d boson 
topological insulator are closely related to those describ- 
ing deconfined criticality of two dimensional bosons. For 
example, a possible direct transition between superfluid 
and valence bond solid (VBS) order controlled by a de- 
confined quantum critical point was discussed iri 2 These 
field theories naturally incorporate the restriction that 
there is no trivial insulating phase. There is however one 
important difference. For the surface states discussed in 
this paper, the trivial insulating phase does not exist even 
in the presence of disorder that breaks translation sym- 
metry. Correspondingly gapless surface critical points 
are described by random versions of deconfined quantum 
criticality. In contrast, the proofs of LSM and its gen- 
eralizations assume translation symmetry since bosons 
on two dimensional lattices can always form a localized 
phase in the presence of sufficiently strong disorder. This 
is consistent with the fact that at least one of the order 
parameters that enter in 2D bosonic deconfined critical 
points involve spatial symmetry (like the VBS order), 
while the unique property of the surface state of a 3D 
topological phase is that it is able to achieve the same by 
invoking purely internal symmetries. 

III. TRANSPORT PROPERTIES OF 3D 
BOSONIC TOPOLOGICAL INSULATORS: 
GENERAL CONSTRAINTS 

We begin by considering a system of interacting bosons 
in d = 3 space dimensions in the presence of time reversal 
and particle number conservation symmetries. Specifi- 
cally let us consider the situation where the boson field 
b carries charge 1 under a global U(l) symmetry and 
transforms as b — > b under time reversal. The corre- 
sponding symmetry group is U(l) X Z% . Assume it is in 
a gapped insulating phase (at least in the absence of any 
boundaries) and that there is a unique ground state on 
topologically non-trivial manifolds. For any such insula- 
tor in 3D, the effective Lagrangian for an external EM 
field obtained by integrating out all the matter fields will 
take the form 

£eff = £-Max + £-6 (3) 
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The first term is the usual Maxwell term and the second 
is the 'theta' term in Eqn. [I] 

Several properties of the theta term are well known. 
First under time reversal, Next on closed mani- 

folds, the integral of -J^E.B is quantized to be an integer 
so that the quantum theory is periodic under 9 — > 9 + 2tt. 
These two facts together imply that time reversal sym- 
metric insulators have 9 — nir with n an integer. Trivial 
time-reversal symmetric insulators have 9 = while free 
fermion topological insulators have 9 = tt. 

If we allow for a boundary to the vacuum and fur- 
ther assume that the boundary is gapped (if necessary 
by breaking time reversal symmetry), then the 9 term 
leads to a surface Hall conductivity of To see this, 
assume a boundary (say at z = 0), 9 = 9(z) is zero for 
z < and constant 9 for z > 0. The action associated 
with the 9 term is 

S e = i J d 3 xdt9(z)d^ (4) 

1 f i d9 , 

= -^j dxdt Tz K ® 

= A / d 2 xdte zvXK A v d x A K (6) 
87r Job 

Where A is the external electromagnetic potential and 
= e^ Xli A v d\A K . This is a surface Chern-Simons 
term and leads to a Hall conductivity 9/2-k. 

For fermion topological insulators 9 = tt so that the 
surface cr xy = \. If we shift 9 — > n + 2nir, then the surface 
&xy = ( n + \)- This corresponds to simply depositing 
an ordinary integer quantum Hall state of fermions at 
the surface of this insulator - hence this should not be 
regarded as a distinct bulk state so that the only non- 
trivial possibility is 9 = 7:. 

Now let us consider bosonic insulators. Again T- 
reversal and periodicity imply 9 = nir and a surface 
o~ X y = n/2. A crucial observation is that now 9 = 2tt 
must be regarded as distinct from 9 = 0. At 9 = 2ir 
the surface a xy = 1. But this cannot be obtained from 
the surface of the 9 = insulator by depositing a ny 2d 
integer quantum Hall state of bosons. Recent worlP^E^ 
has shown (see Ref. [T3]for a simple argument) that 2d 
IQHE states of bosons necessarily have a xy even. Thus 
the surface state of the 9 = 2ix boson insulator is not a 
trivial 2d state but rather requires the presence of the 3d 
bulk. 

Therefore 9 = 2tx necessarily corresponds to a non- 
trivial 3d bosonic TI. 9 = An is however trivial as then the 
surface state can be regarded as a 2d bosonic IQHE state. 
One may still obtain a 3D topological phase, but the 
topology is not manifest in the electromagnetic response. 

We can sharpen and generalize this result. Under T- 
reversal as 9 — » —6, nir — » —nir. As the bulk state is 
T-reversal invariant we require that the surface state at 
9 = —nir be obtainable from the surface state at 9 = +nir 
by depositing a 2d IQHE boson state. Let us charac- 
terize the surface state by both its electrical and ther- 
mal hall conductivities {<r xy , n xy ). Under T-reversal both 



Hall conductivities change sign. The requirement above 
then means that (2a xy , 2n xy ) must correspond to the al- 
lowed electrical/thermal Hall conductivity of a 2d boson 
IQHE state. 

For 9 = 2tt, it follows that a xy = 1, n xy = 0. It is thus 
"half of the elementary 2d boson IQHE state. 

For 9 = 7r, 2a xy = 1 and this is not allowed for the 2d 
bosonic IQHE. It follows therefore that 3d bosonic TIs 
with no "intrinsic topological order" cannot have 9 = ir. 
It is of course very easy to construct such states^^U 
(or other states with fractional 9) if we allow for frac- 
tionalization of the boson but that violates our original 
assumption. 

A 2d IQHE state with 2a xy = 8, 2k xv = 8 is allowed 
and is discussed by Kitaev. Thus a 3d boson TI with 
surface a xy = 4, n xy = 4 is allowed. Combining these 
two types of fundamental states generates the allowed 
thermal and electrical Hall responses on the surface. 

Later on in the paper we will discuss how these re- 
sults fit in with the formal classification of SPT and other 
short ranged entangled phases in 3d. For now we reiter- 
ate the crucial observation of this section: A state with 
EM response of 9 = 2tx necessarily describes a topologi- 
cal insulator of T-reversal symmetric bosons while 9 = it 
requires the presence of "intrinsic topological order". In 
the next Section we study the properties of this 9 = 2tx 
boson topological insulator in detail. 



IV. SURFACE THEORY OF 3D BOSONIC SPT 
PHASES 

In this section we will derive the non-trivial surface 
theory of one example of a 3D bosonic SPT phase. We 
will soon specialize to the symmetries of the topological 
insulator: charge conservation and time reversal sym- 
metry (U(l) xi Z2), and exhibit a nontrivial topologi- 
cal phase in three dimensions, built purely of bosons. 
To begin with we will assume that there are two species 
of bosons whose numbers are separately conserved, and 
there is enlarged (U(l) x U(l)) xi Zj. Later we will 
break this to just U(l) x Zj symmetry by including inter- 
species boson mixing terms in the Hamiltonian. A similar 
const ruction^ has proven to be very powerful in d = 2. 
We consider two approaches: 

The first approach exploits the fact that a bosonic 
SPT phase in d dimensions has surface states that cor- 
respond to a conventional theory of bosons in d — 1 
dimensions except in the way symmetries are imple- 
mented. For example, the edges of SPT phases of 
bosons in D=2+l dimensions correspond to conventional 
ID Luttinger liqui ds, exce pt for their unusual symme- 
try transformations^- 2 * 13 * 37 !. We therefore consider a 2D 
bosonic state to model the surface and assume the sur- 
face is a superfluid breaking one of the U(l) symmetries. 
Then, vortices of this condensate may transform under 
a projective representation of the remaining symmetry 
group. In a projective representation, even the identity 



element of the symmetry group induces a phase rotation. 
Hence local operators, that can be physically measured, 
must remain unchanged under the identity operation of 
the symmetry group, since this corresponds to 'doing 
nothing'. However, vortices, which are non-local objects, 
can transform projectively. One may attempt to restore 
the U(l) symmetry by condensing vortices. However, 
the projective transformation ensures that when vortices 
condense, they necessarily break another symmetry. In 
this way both the boson and vortex condensates lead 
to symmetry breaking, in line with our general expec- 
tation for the surface of a 3D SPT phase. It is important 
that vortices transform projectively, so that they cannot 
be screened by bosons to obtain a trivial representation 
of the symmetry group. This is a generalization of the 
idea of quantum number fractionalization - for example 
a particle with half charge changes sign under the 2ir 
phase rotation of bosons, implying a projective represen- 
tation. Clearly, a half charge cannot be screened by any 
finite number of bosons. Projective representations were 
also recently used to classify SPT phases in D=l+lwhere 
they correspond to the ends of gapped one dimensional 
topological phases. For example, the half-integer spin 
edge states of spin-1 Haldane chains furnishes a projec- 
tive representation of the rotation group. This suggests a 
physical picture of a 3D SPT phase, in which the vortex 
line in the bulk is similar to a Haldane chain type gapped 
phase, which necessitates low energy states on the sur- 
face where the vortex ends. In this section we specialize 
to the symmetries of the topological insulator. Then this 
procedure explicitly produces a topological phase char- 
acterized by quantized magnetoelectric effect 9 = 2ir. 

The second approach will be to directly implement the 



property discussed above in Section III that if the surface 
breaks T-reversal and is gapped then it has quantized 
Hall transport. If T-reversal is not broken a powerful ap- 
proach to obtain the surface theory is to start with the 
theory of the quantum phase transition point between 
the two bosonic quantum Hall phases that correspond to 
the two T-broken surfaces. In the case of free fermion 
topological insulators, a similar reasoning leads to the 
single Dirac cone surface state that describes the tran- 
sition between the a xy = ±| states on the surface. For 
free fermions the transition between these integer quan- 
tum Hall states is described by a Chalker-Coddington 
network modepSI. For the bosonic problem of interest 
here we construct an analogous network model and show 
it leads to a sigma model with a topological term. 

The results of these two approaches are readily seen to 
be connected. In both cases the field theories we obtain 
for the surface have appeared previously in the context 
of deconfined quantum criticality. We discuss the phase 
diagram of the surface states described by these field the- 
ories. When inter-species tunneling is included the vor- 
tices of the two species of bosons get confined to each 
other. The resulting single vortex no longer transforms 
projectively under the physical symmetries. However we 
argue that it is most conveniently viewed as a fermion. 



This precludes the possibility of obtaining a trivial insu- 
lating phase at the surface by condensing vortices. 



A. Surface States and Projective Vortices 

Consider a boson field at the surface with phase de- 
gree of freedom <fii, b\ — e 1 ^ 1 . We assume the bulk is 
insulating and the surface is in the x, y plane. 

The surface theory could spontaneously break a global 
U(l) symmetry of boson number conservation (a surface 
supcrfluid) or stay insulating. More precisely as the bulk 
is always assumed insulating, the vortex line loops have 
proliferated in the bulk. These vortex lines penetrate the 
surface at points, which may be viewed as point vortices 
of the two dimensional surface theory, since there is no 
vortex line tension in the insulating bulk. These point 
vortices are gapped when the surface is a superfluid. If 
instead they are condensed the surface will be in sulat- 
ing. To describe vortices we go to a dual descriptio n^ 35 ! 36 !, 
where we write the density and currents of the boson bi 
on the surface in terms of the field strengths of a gauge 
field ji M = e^ x d u a2x/2n. In particular the density of 
bosons is: rii — (d x a.i y — d y a2x)/^ 7r (the reason for the 
subscript 2 on a will soon be apparent). The boson in- 
sertion operators e 1 ™^ 1 correspond to monopole insertion 
operators, since they insert 2itm magnetic flux. Now, the 
vortices ^2 are particles that couple minimally to the 
gauge field a^- In general there will be multiple vortex 
species that transform into each other under the symme- 
try operation. Let us label them by i, so [^2)1 — "fe- All 
these fields couple minimally to the gauge field. 

Thus we have for the dual surface theory: 



+v^2<) + + 



(7) 
(8) 



where / 2m „ = {d^a-tv - d v a 2p ). 

As argued above, one route to obtaining topological 
surface states is if the surface vortices transform under 
a projective representation of the remaining symmetry. 
Vortices that transform projectively under a global sym- 
metry arc actually not at all unfamiliar. It describes the 
generic situation of two dimensio nal b osons on a lattice, 
say, at some commensurate filling^m These projective 
vortices play a crucial role in the theory of deconfined 
quantum criticality. 

We will only need to consider two component vortex 
fields, ^2 = (^2+, "02-) for a variety of cases consid- 
ered in this paper. The gauge invariant combination 
^la+^2 = ^2+^2- = e^ 2 = b\ then defines another 
bosonic field . Now, Eqn. [7] closely resembles the action 
for a deconfined quantum critical point (no n-compact 
CPi theory with easy plane anisotropy) 2 40 . We will 
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demonstrate how this emerges as the theory for the sur- 
face states, and also, in the next section, describe a three 
dimensional bulk theory that leads to this edge theory. 



Surface States of a Bosonic topological insulator: 



Symmetry U(1)~A Zi 



Condensing single vortices will then break symmetry 
as described below. This is best analyzed by assuming 
separate number conservation of each species of boson, 
7ii, fi2, in which case they can be coupled to an external 
gauge potential A\, A 2 . Then the effective Lagrangian 
EqrJT] reads: 



These are the symmetries of the topological insulator, 
a conserved U(l) charge and Zl^ time reversal symme- 
try. The semi-direct product appears so that the charge 
insertion operator e 1 ^ is invariant under time reversal, 
which involves both cf> — > — <f> and i — > — i. Here we will 
construct the surface theory of a 3D topological phase 
with these symmetries. 

Let us begin with an enlarged symmetry, two species 
of bosons that are separately conserved. Consider a con- 
densate of one species b\ . Vortices in this condensate are 
created by the field ^2- We need to specify the projective 
representation for the vortices ^ 2 and the transforma- 
tion of the bosons b\. The remaining symmetry group 
U(\) xi Z 2 has a single projective representation (PI) 
which acts as follows. Under a U(l) rotation by angle e, 
the fields tp 2 ± ~^ e ±le ^ 2 ip2± and under time reversal: Zj: 
xp2+ ~> ^2- an d "02- — > 1P2+ ■ Or more compactly: 



e' 2 



:U(1) 



*2 -> <T X % 



(9) 



Here the a are the Pauli matrices in the standard repre- 
sentation. Thus the vortices carry charge ±1/2, of bosons 
of the other species. The time reversal symmetry that in- 
terchanges the two vortex fields ensures that the vortex 
charge is fixed exactly at half. It is impossible to 'screen' 
this charge with regular integer charged bosons. 

These transformation laws of course determine how the 
boson operator b\ = e 1 ^ 2 and their density n 2 transforms. 

We also need to specify how the bosons b\ = e 1 ^ 1 and 
density ni transforms. The symmetry transformations 
are: 



(10) 



The conjugate boson numbers therefore transform as 



"1,2 -> "1,2 

m, 2 "1,2 



U(l) 



zl 



(ii) 



A necessary compatibility check is that Eqn. [7] is in- 
variant under the symmetry operation, which can be veri- 
fied for these transformations. For example, time reversal 
symmetry is implemented via Eqn. [9] on the vortex fields 
which is compatible with n 2 (and hence a 2 ) remaining 
invariant while i — > — i under time reversal. Moreover, 
since the bosons carry charge, the monopole insertion 
operators are forbidden. 



a = E 



+ 77— {tnv\d v a 2X ) 2 + -^-A lfl e^xd u a 2 x (12) 

lK\ 111 

Here s = ±. The ip 2 ± a re vortex fields of b\ which carry 
half charge of boson species 2. The flux of the gauge 
field ct\ is precisely the conserved density of species 1; 
hence the last term in the above action where the exter- 
nal probe gauge field A\ couples to this current. This 
action will need to be modified by including all symme- 
try allowed perturbations. We will do this and analyze 
the possible phases below, but as a preview, consider the 
effect of breaking time reversal symmetry by condensing 
just one species of vortex say ip 2 -. Anticipating a single 
charge, consider an external field that couples equally 
to the two charge densities A\ = A 2 = A. The vor- 
tex condensate forces a 2 = \ A, which when substituted 
into the action above yields the electromagnetic response 
-^A^t^xdyAx, which yields a xy — 1 on the surface, indi- 
cating a magneto electric response 9 = 2n, as advertised. 
We now change track and obtain the same surface theory 
from a very different point of view. 



B. Network Model Construction 



The general considerations of Section [ni] showed that a 
time reversal symmetric boson insulator with electromag- 
netic response characterized by 9 = 2ir is in a topological 
insulator phase. This key result relied on the observa- 
tion that if T is broken at the surface to gap it out then 
such a state has a quantized electrical Hall conductiv- 
ity Cxy = ilj an d a thermal Hall conductivity n xy = 0. 
What if T-reversal is not explicitly broken at the surface? 
The surface can then potentially be gapless. What is the 
nature of the resultant theory? To construct this theory 
it is extremely instructive to learn from the example of 
the free fermion T-reversal symmetric topological insula- 
tor. In that case if T is explicitly broken to gap out the 
surface, then we get a xy = When T is unbroken, it 
is possible to get a single massless Dirac cone which is 
exactly the low energy theory of the transition between 
two integer quantum hall plateaus of fermions in d = 2. 
Generically we can tune the chemical potential to move 
away from the Dirac point to get a Fermi surface which 
encloses the Dirac point. 

The familiar free fermion example gives us a crucial 
clue to construct the theory of the T-reversal symmetric 
surface state of the boson topological insulator. First 
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construct the low energy theory of the d — 2 integer 
quantum Hall state of bosons as a potential candidate for 
the gapless surface state of the 3d topological insulator. 
Then add perturbations allowed by symmetry to obtain 
the generic surface theory. 

With this motivation we now study the IQHE plateau 
transition of bosons in d = 2. 

IQHE quantum phase transition of bosons in 
d= 2 

We will study the phase transition using a "network" 
model construction. The idea is to start with the theory 
of the edge state and couple together opposite edges. Let 
us warm up with the familiar example of the IQHE tran- 
sition of fermions from a state with a xy — 1 to one with 
a xy = 0. The model is defined by Fig. [2] It is described 
by the Euclidean action S — J dxdtC with 



12 3 4 



to 



I r ..2 3. 



to 



te" 



c 



with 



(d T - isid x ) • 



,U (ci+iCi + CjCj+i) (13) 



tieven 

si = -i-iy 



(14) 
(15) 
(16) 



The first term is the sum of the actions of a single chiral 
edge mode of the a xy = 1 fermion IQHE state taken 
to propagate in opposite directions for adjacent i. The 
second term describes electron hopping between opposite 
moving edge channels. 

When t e < t a all chiral edge channels are paired with 
partners and a xy — 0. Conversely if t e > t a , then all 
edge channels get paired except the two end ones and we 
get the fermion IQHE state with a xy = 1. The transition 
occurs at t e = t a and is readily seen to yield a single 
massless Dirac fermion in the continuum limit is taken 
in the i direction. 

Let us now repeat this construction for the bosonic 
IQHE transition. The edge theory for the boson IQHE 
state with a xy = 2, n xy = has a pair of counterpropa- 
gating ed ge mod es - one carries the charge and the other 
is neutral^- 2 * 13 * 37 *. It is convenient to write the effective 
action of the edge as an SU(2)\ WZW theory: 



S, 



dxdr^-tr (d^d^g) + iS< 



wzw[g] (17) 



Here g is a 2 x 2 matrix with entries g — 



The 



b\ , &2 are the two physical boson fields which form the 
IQHE quantum Hall state. 

A network model capable of describing the boson 
IQHE state may now be written down and is defined 
by Fig [3] Again we have an array of opposite edge 
channels which are coupled together by boson hopping 

~ Ea=i,2 ( b ta b i+i,a + h.cj cx -tr (gjg i+ i + h.cj . The 



FIG. 2: Network model for the fermion IQHE transition. On 
the top all chiral edge modes are paired to yield an ordinary 
insulator. On the bottom there is an unpaired edge mode to 
yield an integer quantum Hall insulator 



full effective action is then 



S 
So 

s w 



So + Sw + St 

j dXdT 2~X ^ tT i 9 !* 9 * 9 !* 9 * 

i 

i^SiSwzwidi] 

i 

-^2utr (glgi+i + h.cj 



(18) 
(19) 

(20) 

(21) 



with Si and t% as before. If t a 3> t e we get the trivial 
insulator while if t e 3> t a we get the boson IQHE state. 
The transition occurs at t e = t . A low energy theory of 
the transition is obtained by taking the continuum limit 
in the i direction. As the opposite moving edge channels 
have opposite WZW terms they nearly cancel, and it is 
necessary to carefully sum them. Fortunately precisely 
this sum was performed in Ref. I26l where the same model 
arose in a different context. The result is the effective 
D = 2 + 1 dimensional action 



'e// 



d 3 x^-tr (d^d^g) + iir£e[g] (22) 



The second term is a 9 term for the SU(2) matrix valued 
field g in 2 + 1 dimensions corresponding to n 3 [5L/(2)] = 
Z. In the present context our calculation has yielded this 
term at the value 8 — n. 

We of course do not have full SU(2) symmetry rotat- 
ing between bi and 62 in the microscopic system. For the 
time being let us assume that we have U(l) x U(l) sym- 
metry corresponding to separate conservation of the b\ , 62 
bosons. Further let us also assume there is a Z2 symme- 
try interchanging bi and 62- Later we will relax all these 
assumptions. Then the results of Ref. show that the 
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Y 



FIG. 3: Network model for the boson IQHE transition. Each 
edge channel now has a charged chiral mode and a counter- 
propagating neutral mode. The rest is the same as for 
fermions. 



field theory above at 9 — it maps to the self-dual easy 
plane non-compact CP 1 (NCCP 1 ) model. Equivalently 
it also maps onto a model of two species of spacetime 
loops with a phase ir associated with each linking of the 
two loop species. 

The 9 = 7T SU(2) matrix field theory (with the 
U(l) x U(l) anisotropy) or the equivalent easy plane 
NCCP 1 model arise also in the theory of deconfined quan- 
tum criticality in two space dimensions. Remarkably we 
see that the field theories describing the boson IQHE 
plateau transition (and hence the surface states of the 
3d boson topological insulator) are closely related to the 
theory of deconfined quantum criticality. In the previous 
subsection we obtained this connection from a different 
point of view. We will in the rest of the paper explore 
this connection in more detail and generality . For now 
we merely point out that the results of Ref. (see also 
Ref. EHI) show that the 9 = n SU(2) matrix field theory 
in two space dimensions does not have a trivial gapped 
disordered phase. Its phases either break symmetry, are 
gapless or have topological order. This is a hallmark 
of the surface state of a symmetry protected topological 
phase - there is no trivial gapped phase that preserves all 
the symmetries. 

Surface of the Bosonic TIs: Field theories We 

now exploit our intuition about deconfined criticality to 
obtain the theory of possible surface states of the 3d bo- 
son TI starting with the 9 = ir SU{2) matrix field theory. 
We first describe a number of equivalent field theoretic 
descriptions of the surface state paying particular atten- 
tion to the realization of the physical U(l) x symme- 
try. First we note that the SU (2) matrix g is related to 



the physical boson fields bi t 2 through 

% -b* 



fJ 



b\ 



(23) 



Under the global U(l) symmetry, both bosons transform 
with charge 1, i.e 

h, 2 -> b h2 e^ (24) 

We implement time reversal by simply requiring that 

h, 2 -> &1,2 (25) 

In terms of the phases of the bosons, defined through 
&12 ~ e'^ 1 ' 2 , and the conjugate bosons densities 711,2 the 
symmetry transformations are the same as in Eqns. |10|11| 
so that we are indeed describing the same symmetry class 
in the two approaches. We remind the reader that the 
total number ni+n 2 of the two boson species is conserved 
due to the global U(l) symmetry but the relative number 
rii — n 2 is in general not. As promised before we will first 
analyse the theory in a limit where this relative number 
is also conserved (so that there is U(l) x U(l) symmetry) 
and then include interspecies tunneling terms to recover 
the generic case. 
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TABLE I: Symmetry properties of the ipi±,ip2± fields. qi : 2 
are the charges under the two U(l) symmetries associated 
with 61,2 respectively. n„i_2 are the vorticities in the phase 
of 61,2- They can also be viewed as the gauge charge for the 
coupling to the corresponding U(l) gauge fields. The last 
column gives the transformation under time reversal. 

As argued in Ref. 2BJ the 9 term of the SU(2) matrix 
field theory with U(l) x U(l) anisotropy has a simple 
interpretation. It is the phase that is picked up when 
the vortex of the boson 61 is taken around the vortex 
of the boson b 2 - At 9 = ir the two vortices are mutual 
scmions. We may thus readily write down a dual field 
theory in terms of the vortices $i Vi2 v of the two bosons 
&i,2 respectively. This has the structure 

C = C lv + C 2v + Cg + C A (26) 
Ci V = |(d Al - i (a v +/3 v ))$ l „| 2 + (27) 

Ce = — flint uvxdvfox (28) 
L A = iAi^ji^ (29) 



Here [a, is, X, .... represent space-time indices in 2 + 1 di- 
mensions. The cij,i = 1.2 are the usual dual gauge fields 
of the vortex theory. The physical current jip. 01 the 
bosons bi i2 is given as usual by 



Jifj, 



1 

2^ 



(30) 
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We have included external probe gauge fields that 
couple to these currents. The fop are 'statistical' gauge 
fields that serve to impose the mutual statistics of the two 
vortex species through the mutual Chern-Simons term in 
£g. We have also tuned away a chemical potential that 
couples to the total boson number so that the effective 
action is relativistic. We will shortly relax that assump- 
tion. 

In passing we note that recently related models of two 
species of bosons with mutual it statistics have been stud- 
ied numerically through Monte Carlo simulations^!. The 
relevance of these models to the surface of the 3d boson 
topological insulator (and the related 2d boson integer 
quantum Hall transition) should give further impetus for 
such studies. 



C. Synthesis of the Two Approaches 

We now provide a synthesis of the results of the two 
approaches taken in this section. We will rely closely 
on the results of Ref. [26] to provide two alternate field 
theoretic representations of the theory described by Eqn. 



26 Rather than repeat the derivation from Ref. (261 we 
provide a physical description. The tt phase picked up by 
the vortex $ lv goes around the vortex $ 2 « suggests that 
$i„ carries 1/2 charge under the global U(l) symmetry 
associated with species 2 and vice versa (i.e, <J>2„ caries 
1/2 charge under the global U(l) of species 1). However 
the 7r phase is obtained for both charge 1/2 and charge 
— 1/2. We thus should expect that the vortex of either 
species carries fractional charge ±1/2 of the global U(l) 
quantum number of the other species. This expectation 
is formalized by the derivation^. First by doing a duality 
on one species (say 1) we explicitly map to an easy p lan e 
non-compact CP 1 model with action given by Eqn. 12 



If instead we had performed a duality transformation 
on species 2, we would have obtained an equivalent ac- 
tion in terms of the fractionalized fields ipi± related to &i 
through 



(31) 



This action takes the form 

£ 2 = E 



dp - iatip - i^- ) (.-J., 



1 



+ TT~ {tuvxdvctix) 2 + —AzptuvxdvOiw (32) 

with the physical U(l) current of the &2 bosons given 
by jip = ^£pv\d u ai\. Note the obvious similarity of 



Eqn. [32] with Eqn. [12] after interchange of the 1 and 2 
labels. This is a reflection of the self-duality of the easy 
plane NCCP 1 model first pointed out in Ref. l40lThis 
self-duality is obvious when both theories are obtained 
starting with the sigma model or the equivalent dual vor- 
tex theory (Eqn. 26). 



D. Analysis of the Surface Field Theory: Phase 
Diagram, Deconfined Criticality and Fermionic 
Vortices 

Having obtained a field theoretic description of the sur- 
face states we now analyse the phase diagram. The sym- 
metry transformations summarized in Table. [IJenables us 
to deduce the allowed perturbations to the actions above. 
A crucial allowed perturbation are 'chemical potential' 
terms that couple to the boson number \i\n\ + ^2^2- 
Another crucial allowed perturbation is an interspecies 

boson tunneling term —A ^fe{& 2 + h.cj. Let us first dis- 
cuss the phase diagram when these terms are tuned to 
zero. Depending on the question being asked we will find 
it useful to use one or the other of the formulations pro- 
vided above. For clarity of presentation we will however 
use the theory in Eqn. [l2]to the extent that is convenient. 

I. T Breaking States and Quantum Hall Effect 
Consider condensing just one of the vortex fields in Eqn. 

m 



(V> 2 +) £ 0; = 



(33) 



Such a phase clearly breaks T-reversal symmetry, which 
interchanges the vortices. However, it is an insulator 
because the gauge invariant combination ("02+ i ) 2-) = 
(e 1 ^ 2 ) = 0. The transport properties of this phase are 
readily obtained by noticing that the combination + 
is Higgsed. Therefore at long wavelengths we may set 



a 2 



2 



(34) 



Further we may integrate out the field ip2- in Lagrangian 
[Hi] The effective long wavelength Lagrangian for the 
external probe gauge fields then becomes 



-ilT 



AiptpvxduA 



2 A 



(35) 



Defining the 'charge' and 'pseudospin' probe gauge fields 

4 _ A t +A 2 A — A 1 -A 2 

Sir n 1 -*i.s — n 



we get 



{A Cil tp V xd v A c x - A a ^p V \d v A a \) (36) 



This implies that the charge Hall conductivity o xy = — 1 
while the pseudospin Hall conductivity <j s xy — 1. Taken 
together the thermal Hall conductivity n xy = 0. If on the 
other hand we had condensed tp2- without condensing 
■01+ we would have found the time reversed partner with 
a xy = 1 and o% y = -1. 

Consider now adding symmetry allowed perturbations 
to the action. The surface state described above is 
gapped and hence is unaffected by the chemical poten- 
tial terms if they are weak. The interspecies tunneling 
term destroys conservation of pseuodspin (= 111—112) and 
hence a xy is no longer quantized. However the electrical 
and thermal Hall conductivities continue to be well de- 
fined and will have quantized values a xy = ±1, K xy = 0. 



13 



This is exactly what we expected based on the general 
considerations of Section IIIII above. 

Pictorially, these edge states may be understood by 
considering the bulk system on a solid sphere and assum- 
ing that -01+ is condensed on the top hemispherical sur- 
face while is condensed on the bottom hemispherical 
surface. Then along the equator there is a domain wall 
between the two kinds of surface quantum Hall states. 
At this domain wall there will be gapless one dimensional 
states identical to the edge of the 2d boson IQHE state. 
Specifically there is one charged chiral mode correspond- 
ing to the jump Aa xy — 2 across the domain wall and a 
counterpropagating neutral mode which carries the pseu- 
dospin. When interspecies tunneling is added the quan- 
tization of the pseudospin Hall conductivity is not guar- 
anteed but the neutral edge mode is protected so long as 
charge is still conserved. 

II. Superfluid State Let us now consider T-reversal 
symmetric phases. A simple option is 



V ~ V-\ which carries charge 0. Further under time 

reversal 



(37) 



independent of s. This state has < 62 >^ 0. The gauge 
field «2 is Higgsed by the tp2± condensate. Consider for 
a moment the situation where the boson number is in- 
dependently conserved for each species. Then this state 
breaks the global U(l) symmetry associated with b 2 but 
preserves the other global (7(1) associated with b\. We 
will refer to it as SFi. If on the other hand both 02± 
are gapped then they may be integrated out to leave be- 
hind a Maxwell action for a%. Integrating out a\ then 
gives a Higgs mass for A\ so that the global symme- 
try associated with b\ is now broken. We will call this 
SF 2 . These two phases are separated by a phase tran- 
sition that is described by the putative critical point of 
the easy plane NCCP 1 field theory. In general, a chemi- 
cal potential term can also be added which will tune the 
system away from the NCCP 1 critical point. Apart from 
SFi and SF2, we have the possibility of a phase with 
coexistence of the two superfluid orders. 

Inclusion of interspecies tunneling has a more dramatic 
effect. First there is now no real distinction between SFi 
and SF2 phases so that the phase boundary between them 
disappear. More importantly as the relative phase of b\b 2 
can no longer wind, vortices in b\ are bound to vortices in 
62 ■ Note that tp±± are vortices in b 2 and 4>2± are vortices 
in b\. When we bind vortices in b\ to vortices in 62, the 
resulting vortices are created by fields 



(38) 



with s, s' — ±. Note that for s = —s' V ss > carries charge 
under the single remaining global (7(1) while for s = s' it 
carries charge ±1. Thus the vortices no longer carry frac- 
tional charge. V can be obtained as a composite 

of the boson creation operator and the vortex V+_, V-+ 
so that only the latter are 'elementary'. Furthermore 

V-\ can mix with V |_ due to the interspecies tunneling 

term. It follows that there is a unique elementary vortex 



(39) 



Thus in the presence of interspecies mixing there is a 
unique vortex which does not transform projectively un- 
der the global symmetries. Does this invalidate our ear- 
lier analysis? In particular can we now get a trivial insu- 
lator by condensing this vortex? The answer is no. The 
point which we demonstrate below is that the effective 
action for the vortex V in the superfluid phase is not 
the usual one but rather contains an extra Chern-Simons 
term. The presence of this Chern-Simons term has a con- 
venient rough interpretation. It changes the statistics of 
the vortex to a fermion! A simple way to picture this is 
in terms of the description in terms of the vortex fields 



<& v i 2 in Eqn. 26 The 6 term in the sigma model descrip- 
tion means that the two vortices are mutual semions. It 
follows that their bound state is a fermion. 

To put some meat into this picture we start with Eqn. 



26 In the presence of interspecies tunneling the individ- 
ual vortex fields ^\ v ^2v will be confined but a bound 
combination $ CT = <&i v $> 2 v will survive. It is therefore 
necessary to reformulate the action in terms of <& C v It 
is convenient to do so first even in the presence of the 
enlarged U{1) x £7(1) symmetry and later include the in- 
terspecies tunneling. To do this we introduce another 
field $ sl) = $i v &2 V . The resulting Lagrangian takes the 
form 



£ 



£[$ cl ,,a 
% 



■13+] +£[$ s „,a_ +/3_] 



+ —fi+ntuvxdvfi+x - —fi-^vxdvfi-X (40) 



4tt 



47T 



Here a± = a% ± 02, /3± = Pi ± (3%. a± are the usual 
dual gauge fields whose curl gives the charge and pseu- 
dospin current respectively. The most interesting terms 
are the coupling to the gauge fields f3± which have self 
Chern-Simons interactions. Including interspecies tun- 
neling leads to linear confinement of & sv . The effective 
dual Landau-Ginzburg theory of the superfluid then has 
the usual form but with the additional Chern-Simons 
term as promised. 

III. Topological and Other Exotic Orders A wide 
variety of other phases are possible depending on the de- 
tails of the surface interactions. For instance a gapped 
topologically ordered Z 2 liquid is possible and is accessed 
within the present formulation by condensing the paired 
vortex (-02+^2- + h.c.) without condensing any other 
fields. In this situation the full three dimensional sys- 
tem when placed on a solid torus will have a ground 
state degeneracy of 4 coming from the surface topologi- 
cal order. It is interesting to consider the properties of 
this state a bit more and its relationship with the su- 
perfluid state. In the Z 2 topologically ordered insulator 
there is an unpaired vortex ip 2 + ~ V*2- which survives 
as a gapped excitation and which carries physical boson 
charge 1/2 of the £7(1)2 global symmetry associated with 
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the boson b^. We will refer to it as a 2-chargon. Follow- 
ing standard reasoning this phase may equivalently also 
be understood as a paired condensate of Vh+Vh- + n - c - 
Thus there is another gapped excitation corresponding 
to the field ~ ipl_ which, in the present context also 
carries charge 1/2 of the physical boson b\. We will refer 
to this as the 1-chargon. These two chargons are mutual 
semions as expected for Z 2 topological order. Note that 
they have bosonic self-statistics. In the presence of inter- 
species tunneling, a pair of 1-chargons can mix with a 
pair of 2-chargons. Both species of chargons continue 
to exist as independent excitations but now they carry 
charge- 1/2 of the remaining global U(l). Finally the 
bound state of these two kinds of chargons is a fermion 
which does not carry fractional charge. We can take it to 
be charge neutral. It is convenient to regard this neutral 
fermion as the vison, and the two kinds of bosonic char- 
gons as the other two non-trivial quasiparticles expected 
for a Z2 topological ordered state. These transformation 
laws are summarized in TabldlTl 

Now let us consider the relationship to the superfluid 
state discussed above. This will enable us to clarify the 
nature of the vortices of the superfluid state. Coming 
from the superfluid side the Z2 topological state is ob- 
tained by condensing paired vortices. In the presence of 
inter-species tunneling we argued above that there is a 
unique vortex V. The unpaired vortex survives as a fi- 
nite energy vison in this vortex pair condensate. That 
this vison is a fermion ties in nicely with the observation 
that the superfluid vortex V is conveniently regarded as 
a fermion. Thus as the transition to this Z2 insulator 
is approached the vortex statistics becomes well defined 
and becomes fermionic. 

The topologically ordered phase provides a particularly 
simple perspective on why a trivial gapped paramagnet is 
not allowed. Generally to go from a topologically ordered 
insulator to a trivial insulator we must confine the topo- 
logical quasiparticles. For a Z2 gauge theory, this is done 
b condensing one of the three non-trivial kinds of quasi- 
particles (usually dubbed the electric, magnetic and their 
composite) . For the Z 2 topological state that can appear 
at the surface of the SPT phase we are discussing, the 
electric and magnetic particles are both (half)-charged 
under the global U(l) symmetry, and their condensation 
breaks this symmetry. On the other hand the neutral 
topological quasiparticle is a fermion and hence it can- 
not condense. At the same time, time reversal prevents 
one from altering the Chern number associated with this 
gapped fermion. Thus we see clearly that a trivial gapped 
state obtained by confinement from the Z2 topological 
state is not possible at the surface. 

An interesting property of this Z 2 topological ordered 
state is that it realizes symmetry differently from strictly 
2d systems. Such a strictly two dimensional gapped 
abelian insulator may be described within the usual K- 
matrix formulation. For Z2 topological order K = 2a x . 
If both bosonic chargons carry charge 1/2 as we argued 
then the charge vector r = (1, 1). It is then easy to see 



that the resulting topological phase has non-zero elec- 
trical Hall conductivity, and therefore must break time 
reversal invariance. However when realized at the surface 
of the three dimensional insulator a time reversal sym- 
metric Z2 topological phase where both bosonic quasi- 
particles carry charge- 1/2 is allowed. In Appendix |d| 
we collect together the properties of the Z 2 topological 
ordered state for the various SPT phases discussed here, 
and show using the results of Ref. 24 that they all realize 
symmetry differently from what is allowed in strictly 2d 
systems. 

While such interesting topologically ordered (or other 
even more exotic) states are allowed they are not re- 
quired: the surface could be in a superfluid or T-broken 
insulating quantum Hall state with no ground state de- 
generacy. The most important conclusion however is that 
a trivial gapped insulating state which preserves all the 
symmetries and has no topological order is not possible 
on the surface. This is a key property of a symmetry pro- 
tected topological phase and is satisfied by our example. 



Excitation 


Charge 


r 2 


Boson 1 (e) 


q=l/2 


+1 


Boson 2 (m) 


q=l/2 


+1 


Fermion (f) 


q=0 


+ 1 



TABLE II: Symmetry properties of the topological excitations 
of a Z2 gauge theory, realized on the surface of a 3D SPT 
phase with charge conservation and time reversal symmetry 
(U(l) xi Z% ), as for the topological insulator. The bosonic 
quasiparticles carry half charge, but transform linearly under 
time reversal. Further details are in Appendix |D| 

Some remarks are now in order. We have constructed 
one topological phase, which when the surface sponta- 
neously breaks time reversal symmetry leads to a quan- 
tized magneto electric effect. We do not attempt here to 
construct all possible topological phases in this symmetry 
class, which, according to the cohomology classification 
of Chen et aP gives Zi. Of the three non-trivial states, 
one must have vanishing magneto-electric effect since this 
is an additive quantity. Explicit models for such phases 
is left to future work. In this particular example, the fact 
that there is a single U(l) charge conservation symmetry 
introduces additional terms that imply a deformation of 
the deconfined criticality action. Despite this the degrees 
of freedom of this action provide the useful fields in terms 
of which the effective theory of the surface state may be 
described. We will discuss other examples with different 
symmetries in subsequent sections. 



V. 3D TOPOLOGICAL FIELD THEORIES 

We would like to write down a 3D theory of particles 
and their vortices. We can choose to represent parti- 
cle currents by — e^ Xr7 d u B\ a /2ir. The vortex lines, 
being loops in three dimensional space, sweep out a sur- 



15 



face in space time denned by the two form 7^. Relat- 
ing this to a vector potential a, whose curl is the loca- 
tion of the vortex loop, we define — e fl " Xa 'd\a ( j/2ir. 
The quantization of boson particle number to integers 
implies that the charges that couple to this vector po- 
tential a are quantized. (Equivalently, the dual vec- 
tor potential is compact, i.e. only defined modulo 2tt). 



Clearly, gauge transformations a M 



B^ v — > B^ + {d^a v — d v OLfj) do not change physical 
variables. Taking a particle around a vortex leads to a 
phase of 27r, which is captured by the minimal coupling 
C = a^j^ , which may be rewritten as: 



£ 



1 

2^ 



a^ vXa d v B Xa 



(41) 



This is often written as eBda/2ir, and called the BF 
action. The unit coefficient in the action ensures that 
there is no topological order, i.e. a unique ground state 
in the absence of surfaces^ - appropriate to the current 
discussion. Then, the theory above only states the obvi- 
ous, that particles and their vortices have a mutual phase 
factor of 2ir. 

Let us briefly describe how this term arises from a mi- 
croscopic theory of bosons. A lattice regularized theory 
of bosons can be captured by a loop model of integer val- 
ued closed loops with Euclidean Lagrangian Ct = ^J/Jm- 
The integer constraint is implemented by summing over 
the auxiliary vector field that is an integer multiple 
of 2tt. Now, the current j M takes real values, and its di- 
vergence free condition can be implemented by writing 
iii = ^^u\adijB\ a /2tt , where the two form B is also a 
real field. This gives: 



f E 



1 



8lT 2 



-{edBf 



2tt 



-eadB 



(42) 



where the second term is the desired statistical interac- 
tion (the factor of i appears because of the Euclidean 
formulation). However, at this point a is an integer 
(times 2ir) field. One can softly introduce this constraint 
by assuming a to be real, but adding the cosine term 
AC = — Acos(9 M — a M ), where we utilized the fact that 
longitudinal component can always be added to the gauge 
fielcP^. The phase <fi is actually the phase of the origi- 
nal bosons, and when the bulk is insulating the cosine 
is irrelevant, since the bosons are gapped. Therefore in 



the insulating phase we may use Eqn. 42 where both 
B, a fields taken real, with the caveat that charges are 
ultimately quantized. 

As shown in appendix A, surface states defined from 
the BF theory Eqn. [42] are usual 2D bosonic modes that 
are not topologically protected. To encode a SRE topo- 
logical phase an additional term must be added as shown 
below. 

Based on the discussion on surface states, where a pair 
of bosonic fields were invoked, we consider two species 
of bosons to write down a topological term. This also 
follows the two component U(l)xU(l) symmetric Chern 



Simons approach for the 2D systems, which was found to 
be successful in describing 2D SRE topological phases^. 
Therefore we will introduce two B fields that represent 
their conserved currents, and two a fields which are vor- 
tices in these fields. Let us begin with the general case 
of N species of bosons, with: 



1 N 

Cbf = ^Y, eBidai 



(43) 



i=i 



where e is the antisymmetric symbol and indices have 
been suppressed. 

Note, the apparently more general version is Cb = 
Hj^-eBidaj with Q a uni-modular i.e. detQ = 1 inte- 
ger matrix which ensures absence of topological order. 
Ho weve r, this can be brought into the canonical form of 

[Q- 



Eqr 
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by redefining Bj = [Q~ ]kiB' k . The transforma- 
tion matrix Q~ Y is also an integer matrix, since det Q = 1 
and the minors of an integer matrix are also integers. 
Now, an additional topological term can be added to 



the action: 



C 



3D 



Cff = 



e 



Cff 
Kuedaidaj 



(44) 



The action must be invariant under + 27T, to allow 
for addition of 2D layers at the surface. Therefore, the 
action defined on a closed three dimensional space should 
be invariant under the shift — > + 27r. It is shown in 
Appendix [B] that this condition fixes the entries K]j to 
be integers. 

A stronger condition on K can be applied as follows. 
Values of that differ by 2ir simply correspond to differ- 
ent ways of terminating the surface. Hence, at a domain 
wall where changes by 2tt at the surface, we demand 
that all excitations present are bosonic. This is the same 
as the requirement placed on K matrices describing 2D 
SPT phases, i.e. that DetK — 1 and all diagonal entries 
are even integers. 

The simplest choice of K matrix with these properties 



K = 




(45) 



We find that for most of the 3D bosonic SPT phases 
that we will be interested in, it will suffice to consider 
this matrix. This is similar to the 2D situation^, where 
the above K matrix describes a large set of SPT phases, 
which differ in the way symmetry is implemented. 



A. Two component BF theory of Bosonic 
Topological Insulator 

Let us specialize to the symmetries of the topological 
insulator U(l) xi Zj, and consider a two component the- 
ory with the simplest allowed K matrix given by Eqn. 
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1451 Then we can write: 



C-tot 


— £bf + £ff - 


Cera 


(46) 


£>bf 


= —e(Bxdax-\ 

ITT 


- B 2 da 2 ) 


(47) 


Cff 


Att z 




(48) 
(49) 



We will discuss coupling to the external electromagnetic 
field C em subsequently. Under time reversal symmetry 
we have: 



(50) 
(51) 



where indices i, j refer to spatial coordinates. The trans- 
formation of B fields is obtained by relating them to the 
boson densities and currents, while the a fields are cho- 
sen to transform such that the BF term is left invariant. 
Since both species aj transform in the same way under 
time reversal we may conclude, — > — under Z 2 . 
A time reversal invariant bulk action can then be con- 
structed for = 0, tt, (given the ambiguity in modulo 
2tt). We of course pick = tt in the topological phase. 
While we have not derived this action, we have writ- 
ten down the simplest possible topological theory which 
meets the general constraints required of SPT phases. We 
now proceed to show it produces a surface with the same 
physical properties as predicted in the previous section. 
We study three different situations: first, we study the 
surface superfluid, and determine the quantum numbers 
of vortices. Second, we investigate the electromagnetic 
response, particularly the magneto-electric polarizability. 
Finally, we analyze the case where time reversal is bro- 
ken at the surface, in opposite ways, leading to a domain 
wall. 

(i) Fractionally Charged Vortices: Consider a surface 
of the topological phase at z = with = tt (0) for 
z > (z < 0). Then the effective action at the surface 
arising from Cff is: 



1 

47T 



dtdxdy e zal31 a\ a dpa 2l 



(52) 



where indices a, ft, 7 run over t,x,y, and the fields are 
evaluated at z = OPP. At the surface we may replace 
a>n — di<pi (see appendix [Aj, where i — x, y. Consider 
now a surface superfluid of component 1 — 2, with a vor- 
tex at the origin x = y = 0. This implies <f> 2 winds around 
the origin, or that (d x d y — d y d x )(f> 2 — 2irn 2 5(x)5(y), 
where we have allowed for a vortex of strength n\. Sub- 
stituting this in Eqn. [52] we have: 



S a = ^ / <///>2<V>i 



(53) 



similarly, a vortex in the field (fix couples to the phase 
of <j) 2 . Given the conjugate relation between number and 



phase, this implies that a vortex of strength n 2 in com- 
ponent 1 — 2 carries charge n 2 /2 of component 1 = 1. 
Thus we see that unit vortices in one bosonic field carry 
a half charge of the other field. 

A different perspective on this result is obtained by 
thinking about the fate of 'external' monopoles of the 
gauge fields a\,a 2 . These are sources for vortex lines 
of the two boson fields bx and b 2 respectively. Now the 
well known Witten effect implied by the term tells us 
that a 2tt monopole in a\ carries gauge charge ^- = g 
that couples to a 2 , i.e it carries charge 1/2 of the global 
U(l) associated with b 2 . Similarly a 2tt monopole in a 2 
carries a charge 1/2 of the global U(l) associated with 
&i< At the surface this monopole creates a 2ir vortex of 
the corresponding boson which is then seen to carry half- 
charge of the other boson thereby recovering the result 
above. 

(ii) Electromagnetic Response: Since one has U(l) 
symmetry, we can couple to an external electromagnetic 
field and write the following terms: 



c = A 



("fidai + j 2 da 2 )edA 



2tt 



{q x Bx +q 2 B 2 )edA 



(54) 



Note, in the second pair of terms qi is the electromag- 
netic charge on the 7th boson given that the current of 
bosons edB = j couples minimally to A. The first term 
is allowed by symmetry and, we will see, required by 
the fact that surface vortices are charged. Since a and 
A transform in the same way under time reversal is 
again quantized to 0, tt, and we will need the latter for 
consistency. 

One can relate qi and 7/ by the following argument. 
At a surface where changes by (50 one can see that vor- 
tices in the two boson species ending on this surface carry 



charge |S(7i, 72)- However, we just noted from Eqn 53 
that the vortices also carry one half of the quantum num- 
ber of the other boson species, which are parameterized 
by Qi, 92- Therefore, when SO — tt we demand 



7i = 92, 72 



(55) 



More formally, we may identify the boson currents via 
j^- [Cbf + Cff], and then coupling to the exter- 
qijjA^, which yields the same result as 



Eqns. 54 and 55 



Lem " 4^ 2 



(q 2 dai + q\da 2 )edA 



2tt 



(qxBx +q 2 B 2 )edA 



(56) 



To find the electromagnetic response, we integrate out 
the Bs, and then the a. The first step gives aj = —qiA. 
Substituting this we get: 
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C 



e 



em— response 



4tt 2 



qiq2\e.dAdA 



Setting 6 = w,q\ = q-x = 1 this is £ em 



; - 



(57) 
(58) 

E ■ B. 



This corresponds to a magnetoelectric polarizability 9 = 
—2ir = 2ir (mod47r), i.e. it is an odd multiple of 2ir, as 
expected. 

Surface Domain Wall: Finally, we consider an insu- 
lating surface in the z — plane, on which time rever- 
sal is broken in opposite ways for y > 0, y < 0. This 
leads to edge modes along the x direction, localized near 
z = y = 0. This is modeled with a spatially varying 
field, where Q(z > 0) = 0, while 6(z < 0, y) = 7rsign(y). 
Introducing this profile in Ctot, with gapped B field on 
the surface, the edge theory is readily shown to be: 



S, 



domain— wall 



dtdx [d x 4>\ <9 t 2 



2tt 



+A d x (q 2 (t>i + qifa) + 



(59) 
(60) 



where the first term defines the commutation relations of 
a regular Luttinger liquid, the second term identifies the 
coupling to the external field (assuming a gauge where 
A x = 0), and the dots refer to non universal potential 
terms for the edge fields. This is identical to the edge 
state of the Integer quantum Hall effect of bosons^EEl, 
on setting qi = q 2 = 1, which is also constant with the 
magneto electric polarizability of 9 = n in this phase. 

The general problem of deriving the 3D field theory 
above, from microscopic models is left for future work. 
Below, we describe a bulk non-linear sigma model which 
is some ways may be considered a microscopic theory 
since it assume additional ingredients, over and above 
the purely symmetry group based approach of the co- 
homology theory^. Therefore, instead of writing sigma 
models with topological terms, where the target manifold 
is the symmetry group, we will allow the target mani- 
fold to be the four sphere S 4 , which assumes a partic- 
ular microscopic representation. However, since we are 
not concerned with classifying phases, but rather with 
providing physical examples, this additional assumptions 
will be convenient. This is quite analogous to the com- 
mon practice of considering crystalline band structures 
for free fcrmion topological insulators, although they can 
(and strictly speaking should) be defined in the absence 
of translation symmetry^. 



B. Bulk sigma model theory 

The bulk field theory discussed in the previous section 
is topological and has no bulk dynamical degrees of free- 
dom. In this section we describe a different bulk theory 
with dynamical boson fields which gives rise to the topo- 
logical effective field theory of the previous section in a 



disordered phase. This theory may thus be viewed as a 
field theory realization of a model with a bosonic SPT 
phase in three space dimensions. This field theory takes 
the form of a 3 + 1 dimensional non-linear sigma model 
supplemented with a topological 9 term. This general- 
izes to three space dimensions the continuum field theory 
model that realizes the 2d integer quantum Hall state of 
bosons. 

Following the 2d example and the discussion in previ- 
ous sections of this paper we will enlarge the symmetry of 
the boson system from U(l) to a larger symmetry group 
and then add perturbations to reduce to the symmetry 
of interest. For the construction of this section it is ex- 
tremely convenient to consider a generalization where we 
first embed the U(l) symmetry into an 5*0(5) group. 
Consider therefore a 5-component unit vector field n. 
Later we will describe exactly how the physical symmetry 
(U(l) x or U(l) x Zj , etc) are realized by the compo- 
nents of this field. For now we write down a continuum 
field theory for h. On a closed space-time manifold, say 
the four sphere, the Lagrangian takes the form 



C 

C 6 [h] 



Co[h] + £s[n] 
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(61) 
(62) 
(63) 



Q = — J d 3 xdTe abcde n a d x n b d y n c d z n d d T n e (64:) 
= — - j d 3 xdr det{h d x h d y h d z h d T h] (65) 



o 2 

where ^4 = is the volume of the unit four dimen- 
sional sphere. Q is the integer invariant corresponding to 
n^S 4 ) = Z and counts the number of times spacetime 
configurations of the n field wrap around unit 4-sphere. 
Clearly the theta term does not affect the bulk physics if 
9 = 2mr with integer n. We note that through out this 
section 9 will denote the theta angle for the bulk sigma 
model, and should not be confused with the same symbol 
used previously for the electromagnetic response. We are 
interested in disordered phases of this field theory where 
n is gapped in the bulk. 

Consider a spatial domain wall between a state where 
9 = 2ir and one where 9 = 0. If this domain wall has a 
non-trivial surface state then the 9 = 2ir theory describes 
an SPT phase in the bulk. Such a domain wall corre- 
sponds to a situation where 9 varies spatially. To handle 
this it is convenient to elevate 9 to be a new dynamical 
field and define a sigma model for a new 6-component 
unit vector field 4> defined by 



cos a 
n sin a 



(66) 



The field <j) defines a map from spacetime (taken to be the 
four dimensional sphere S* 4 ) to the five dimensional unit 
sphere S 5 . Consider a field theory for 4> which includes 
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(apart from the usual gradient terms) a Wess-Zumin- 
Witten (WZW) term denned as usual as 2ir times the 
fraction of the volume of 5 5 that is bounded by the hyper- 
surface traced out by <j>. Formally let <j>(x, u) be a smooth 
extension of <p(x) such that (p{x, 0) = 4>q, <p(x,l) = <f>. 
Then 



Swzw[4>] = ~2 / / dudet[4> d x <p d y <p d z cp d T <j> d u (p] 

n Jx,T JO 



Consider an extension where 



cos a(u) 
n sin a(u) 



(67) 



(68) 



with a(0) = 0, a(l) = a, a(u) independent of x, r) and h 
is independent of u. The determinant in Swzw is readily 
calculated and reduces to the theta term Cg [n] for the h 
field with 



10 



da' f 



(69) 



In particular when a = we get = and when a = tt 
we get = 2-7T. Thus the WZW model for with constant 
a = 7r describes the 3 + 1-d non-linear sigma model for 
n dX 9 — 2ir. 

To study a domain wall between 9 = 2ir and = 
along, say, z = let a — a(z) such that 



a(z — > — oo ) = 

a(z — > OO) = 7T 



(70) 
(71) 



Further we assume that 5 s is localized to within a short 
distance of z = 0. To evaluate the WZW term in this 
configuration it is convenient to use a different extension 
of the 4> field. Specifically let 



cos a(z) 
h(x, r, u) sina(z) 



(72) 



with a now independent of u,x,y,T and n(a;, 0) = no, 
n x ,\ — n(x). The determinant in Swzw is again readily 
evaluated and becomes 



5< 



WZW 



dz 



da 

dz 



det[h d x h d y h d z h d T n d u n] 
(73) 

As is localized at the domain wall at z — 0, we can 
replace h in the integral by its configuration at z = 0. 
The z-integral can now be performed and leads to 



C 



WZW 



Air 



det[n d x fi d y n d z n d T h d u n] (74) 



This is exactly the WZW term (at level 1) for the n field 
at the boundary. Thus the domain wall in question is 
described by a 2 + 1 dimensional 5*0(5) non-linear sigma 
model with a WZW term. 



This field theoretic result is very useful to construct 
a bulk sigma model description of the SPT phases dis- 
cussed in this paper. The simplest application is to 
bosons with symmetry U(l) x discussed in detail 
in the next section. To illustrate this let us first in- 
troduce an U(l) x SO(3) anisotropy and write n — 
[Reip, Imip, N], where under the global U(l) symmetry 
choose ip —> e ie ip, but N —> N. Under time reversal, we 
let V -> i>*N -> -N. Finally under the global SO(3) 
symmetry N transforms as a vector while ip is invariant. 

The level- 1 WZW term plays the following crucial 
role^ in this field theory: it implies that the vortex of 
the ip field transforms as spin- 1/2 under the SO (3) rota- 
tion. Indeed the 5-component sigma model with global 
50(3) x U(l) and symmetries implemented this way, 
and supplemented with a level- 1 WZW term precisely 
arises also as the theory of the deconfined critical point 
between Neel and VBS states in 2d. There the spin-1/2 
attached to the vortex captures the physical picture that 
a VBS vortex is a spinon. 

For our present purposes we need to further explicitly 
break the 50(3) symmetry while preserving time rever- 
sal. Then the N field is no longer a freely fluctuating 
variable. However the crucial point is that as the vor- 
tices of ip form a spinor their Kramers degeneracy will 
be preserved so long as is preserved even if the full 
50(3) is not present. This is exactly the defining prop- 
erty of the surface theory of one of the SPT phases for 
bosons with i7(l) x Z^ symmetry described in the next 
section. We have thus obtained this surface theory from 
a bulk sigma model. 

(i) Meaning of theta term of bulk sigma model 

The sigma model description is useful as it suggests 
a route to obtaining a physical realization of this SPT 
phase. First let us understand the meaning of the bulk 
9 term in this sigma model. As a topological term it 
depends on the global configuration of the n field. In 
general for a theory of TV, ip with 50(3) x U (1) symmetry, 
if there are no topological defects in either the N or the ip 
field it is easy to see that the 9 term vanishes. In 3d the 
TV field admits point hedgehog defects while the ip field 
admits vortex loops. The 9 term implies that during 
a process where a hedgehog is taken around a vortex 
line a phase e accumulates. For 9 = 2ir it follows that 
the hedgehog of N has charge 1 under the global U(l) 
symmetry associated with ip (this ensures that it acquires 
phase 2tt when it moves around a vortex line). Thus 
this kind of SPT phase may potentially be engineered 
by constructing a physical situation where the hedgehog 
defect of a 3-component order parameter is charged under 
the global U(l) symmetry of the bosons of interest. 

(ii) Other SPT phases 

To describe the other SPT phase discussed above with 
the same symmetries, we must implement the symmetry 
differently. First we consider a different situation where 
the 50(3) vector only has U(l) x Z^ symmetry; then we 
write h — (Reipi, Imipi, Reip2, Irmp2, N z ). Following the 
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previous sections we take both ip\ and ip2 to be charged 
under the global £7(1) symmetry. For £7(1) x we take 
under time reversal V>i,2 ~ ^ "01 2 an( i N z — > —N z . If time 
reversal is preserved, then (N z ) = 0. In the presence 
of anisotropy that favors the n to have zero component 
of N z we may drop it to obtain an effective field the- 
ory for the surface. Then the level- 1 WZW term for h 
becomes the familiar 9 term at 9 — n for the remaining 
4-components in the 2 + 1 dimensional surface theory. To 
understand the bulk, note that when there is £7(1) x Z^ 
anisotropy on an 5*0(3) vector the vortex lines of the 
£7(1) field ip2 come in 2 kinds which are distinguished by 
the sign of N z in the core. The hedgehogs of the orig- 
inal N field are then domain walls within the cores of 
these vortices where the N z changes sigrP^l. Thus the 
hedgehog must be regarded as a composite of two kinds 
of monopole sources for the two kinds of vortex lines. 
Formally we may write the hedgehog creation operator 
fv as 

h) = roL m2- (75) 

where m 2 ± create the two kinds of monopole sources. 
Now the charge 1 of the hedgehog implied by the bulk 
9 term implies that these monopole sources are charged. 
Further as Z% changes the sign of N z , it interchanges the 
two monopole sources. It follows that the monopoles m| 
each carry charge ±1/2 of the i/ji field. This is precisely 
what is implied by the Witten effect as applied to the 
two component BF + FF topological field theory of the 
previous section. 

Exactly the same description can also be provided for 
the boson topological insulator with symmetry £7(1) x 
Z 2 ■ Then we consider anisotropy similar to the above 
with tpi : 2 charged under the global £7(1) symmetry but 
let ipi,2 -» ipi,2, N z — > — N z under Zj. The rest of the 
discussion is identical to the one above. 

This establishes the connection between the bulk sigma 
model and topological field theory descriptions. Apart 
from giving an alternate perspective we hope that the 
ideas of this section will provide insights on physical re- 
alization of these SPT phases, a task we leave for the 
future. 



VI. OTHER SYMMETRIES: 

Now let us study various other symmetries. In partic- 
ular we highlight two cases: (i) ZT± time reversal sym- 
metry. This is the simplest symmetry that produces a 
topological phase and we construct a nontrivial phase 
thus indicating a Z2 class. In the absence of a conserved 
charge there is no quantized magnetoelectric effect. A 
separate topological phase with this symmetry, but chi- 
ral modes on a domain wall is constructed in the next 
section. (ii)£7(l) x Z^ This corresponds physically to 
a time reversal invariant spin system where the z com- 
ponent of spin is conserved. Two nontrivial phases are 



constructed. The first has a quantized 9 = 2ir, but the 
symmetry prohibits background charge which allows us 
to sharply define statistics of vortices. The possibility 
of an exotic type of bose liquid, the vortex metal, as a 
generic surface state is discussed. The other nontrivial 
phase has 9 = 0. However, in this case we show that a de- 
confined quantum critical action could emerge on tuning 
just a few parameters. We relegate to the Appendix the 
following symmetry which is also readily analyzed: (iii) 
£7(1) x Z2 for which we obtain Z2 topological phases. 
This is of interest since it does not involve time reversal 
symmetry. 

A. Symmetry Z% 

We now consider the case of just time reversal symme- 
try, both by analyzing the projective representations of 
surface vortices and by constructing bulk field theories 

Surface Theory: As usual it is convenient to as- 
sume a slightly bigger symmetry to identify the relevant 
physics and then break it down to the physical symmetry. 
Here it is sufficient to enlarge the symmetry to £7(1) x Z^ 1 
so we may discuss vortices in a boson field b\ = e 1 ^ 1 . 

Let us first discuss the transformation of the <j)\ field 
under time reversal. If this was like 'charge-phase', then 
under time reversal <j)\ — » —<f>i + r/n, where 77 = 0, 1. 
However, we would be able to pin this for either value of 
r\ by either a cos (f>i or sin <\>\ term. So this does not corre- 
spond to an SPT phase boundary state, since the surface 
can be gapped in a trivial fashion, without breaking sym- 
metry. The other option is that cf>i transforms like the 
XY spin, i.e. (f>\ — > 4>\ + rjir. In this case, for 77 = 1, the 
term that can be added to the Lagrangian is cos(20i +c) 
and <f>i field cannot be gapped without breaking symme- 
try. Now we need to consider the transformation of the 
vortices because that can lead to a gapped state even if 
4>i itself cannot be condensed. Before doing that we note 
that the field rii conjugate to the phase (f>i transforms 
as ni — > —ni under time reversal (i.e. like S z spin) to 
preserve the commutation relations. 

Let us now discuss the transformation properties of 
the vortices, under the remaining symmetry Z^. What 
are the projective representations of the symmetry group 
Zj? These are essentially the end states of a ID topo- 
logical phase with this symmetry. It is well known that 
there is a Z 2 classification of such phases and the non- 
trivial phase is just the Haldane phase with gapless edge 
states which are spin 1/2 objects. Therefore the projec- 
tive representation of Z^ which we need is the following 
transformation of two vortex fields (just like spin 1/2) 
under time reversal, so ip2+ ~^ +1P2- , "02- ~"02+ or 
more compactly 

* 2 -> io y ^2 (76) 

Note, that if the vortex fields condense, they can either 
condense individually, or simultaneously. In the former 
case the gauge invariant operator N z = |'02+| 2 — 1^2- 1 2 
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takes on a nonzero expectation value. This object, under 
time reversal, transforms as N z — > —N z , so this breaks 
time reversal. If both condense, then the gauge invariant 
field 1^2+^2- = e 2 * 2 acquires an expectation value. How- 
ever, under time reversal this also transforms nontrivially 
02 02 + 7T (note and 2 transform the same way) 
and cannot take on an expectation value without break- 
ing time reversal symmetry (both these are spin opera- 
tors essentially). A third option is that ip2s do not con- 
dense individually, but a pair condenses. This leads to a 
Z 2 topologically ordered state that does not break sym- 
metry, which is also consistent with a topological surface 
state. 

The effective theory for bosons at the edge consistent 
with this symmetry is 

£ e = ^ |(d M - ia2^)ip2s\ 2 

S 

m 

where V%rn = e 12 " 1 ^ 1 is the 2m monopole insertion opera- 
tor, which is allowed once we break the U(l) symmetry to 
leave just the time reversal invariance which allows even 
numbers of monopoles. Note, the background magnetic 
field (d x a2 V ~ d y a2 X ) — 27TOi, is odd under time reversal 
symmetry and is not allowed. Note however the boson 
mixing terms e 1 ^ 1 ^ 21 are allowed by symmetry. Other 
symmetry allowed terms are the same as with U(l) x 
symmetry which is discussed below Eqr ]81| 

Surface Z 2 Topological Order and Symmetry: Just as 
we did in Section [IV D| it is extremely instructive to con- 
sider the question of why a trivial paramagnetic state 
is not allowed at the surface from the point of view 
of the Z 2 topological surface state. Here of the three 
non-trivial topological quasiparticles, the two bosonic 
ones each transform as Kramers doublets under Zj , i.e. 
T 2 = — 1. They are simply the unpaired vortex of either 
&i or 6 2 . They have mutual semionic statistics. Their 
bound state is a singlet (T 2 = +1) but it is a fermion. 
To destroy the topological order we must condense one 
of these non-trivial quasiparticles. However when cither 
of the two bosonic excitations condense Z^ is sponta- 
neously broken. The fermion cannot condense, and time 
reversal symmetry prohibits nontrivial Chern number for 
fermions. Thus there is no possibility of a trivial param- 
agnet. 

The cohomology classification^ also produces Z 2 
phases with this symmetry. Our analysis gives a direct 
understanding of the allowed surface structure of the one 
non-trivial member of this symmetry class. 

3D Bulk Lagrangian: Consider the bulk La- 
grangian: 

C BF = ~(B 1 da 1 +B 2 da 2 ) + G~da 1 da 2 (78) 

the first term is invariant under time reversal if we as- 
sume Bf -> Bf and a] -> a} while B\ j -> -Bf and 



oPj — > —a® under time reversal. The transformation law 
for B is obtained by assuming it is connected to a con- 
served current that transforms like spin current. Thus, 
unlike Chern-Simons in 2D, the bulk action is naturally 
invariant under T. However the second term changes sign, 
if both ai transform in the same way. This fixes O to one 
of two values and 7r yielding at least two phases. Since 
there is no conserved charge there is no coupling to an 
external field. However, the 'fractionalized' degrees of 
freedom at the ends of vortices is captured in this for- 
malism. 

It is relevant to note that phenomena used to dis- 
tinguish topological phases previously do not apply for 
this symmetry. For example, the absence of a conserved 
charge does not allow us to define the magneto electric 
polarizability. Also, it turns out that domain wall be- 
tween opposite time reversal symmetry breaking surfaces 
do not carry gapless modes. Recall, in free fermion 3D 
topological insulators and class Dili topological super- 
conductors there are chiral edge modes on T-breaking 
domain walls. Also, for bosonic topological insulators 
there was a non chiral but protected domain wall mode. 
Here the domain wall is nonchiral (as can be seen from 
the K = cr x matrix that enters in the second term) and 
in the absence of a conserved charge the oppositely prop- 
agating modes can acquire a gap. Nevertheless, the sur- 
face states are still special, and are either gapless, break 
symmetry or develop topological order. 



B. Symmetry U(l)x Z% 

Here the U(l) can be interpreted as spin rotation 
(rather than charge conservation) symmetry about z axis. 
However, we will continue referring to the conserved 
quantity as charge. We construct two different topolog- 
ical phases (and the composition of these phases defines 
a third non-trivial phase), which are interesting for dif- 
ferent reasons. 

Surface Theory of Phase 1 and Deconflned Crit- 
icality action: Here we do not enlarge the symmetry. 
Let bosons b\ = e 1 ^ 1 be charged under the U[l) symme- 
try so: 

0i ->■ 0i +e : 17(1) 
m ni : U(l) 

4>\ -> 0i + n '■ z I 

m -> -ni : Zl (79) 

Next, we consider the vortices ^2 of the field e 1 ^ 1 , and 
specify their transform under the remaining time reversal 
symmetry, which has a single projective representation: 

* 2 -> icr„* 2 : Zl (80) 

. Now, since V^+V^- ~ e 1 ^ 2 we have 2 — > 02 + tt under 
time reversal. The effective field theory is written as: 
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£ e = y2\( d ^- i0l 2^2a\ 2 + ^-f^u 
a 

-A[(^ + ^ 2 -) 2 + h.c.]+V(|vI/ 2 | 2 )) (81) 

The second last term is cos 202 which is preserves time 
reversal symmetry. The flux (d x a2 y — d y a2 X )/'^' n = n i 
vanishes on average since the density n\ changes sign un- 
der time reversal. No monopole insertion operators are 
allowed since changing the flux corresponds to inserting 
conserved U(l) charge. We note that this action is very 
similar to the easy plane non compact CPi (NCCPi) ac- 
tion, proposed as the critical theory between a spin 1/2 
easy plane Neel antiferromagnet and Valence Bond Solid 
(VBS) order. The flux here is just the spin density, while 
the vortex bilinear e 1 ^ 2 correspond to the VBS order. 
In contrast to the square lattice with four fold rotation 
symmetry, here the square of the VBS order parameter 
is allowed, as on a rectangular lattice. An important dis- 
tinction from previously discussed deconfined criticality 
is that here translation symmetry is not invoked. 

Symmetry actually permits other terms in this action, 
for example, the linear derivative terms (dta2i — dia2o), 
which corresponds to electric fields (spin currents) in the 
ground state. Similarly, finite gauge charge is also al- 
lowed in the ground state, corresponding to finite vortex 
density, since vortices here do not break time reversal 
symmetry. This introduces linear time derivative terms 
in the action above. However, if we expand the symme- 
try to include a Z 2 that reverses the orientation of the 
U(l) rotation axis, i.e. that sends n\ — > — n\ (which is 
a rotation by 7r around the S x axis in spin notation), 
then this prohibits the additional terms discussed here, 
since both electric field and gauge charge are odd un- 
der this 7*2- Thus, for the topological surface state with 
symmetry [U(l) xi Z 2 ] x Z 2 the field theory is given by 
[8~T| Parenthetically we note that precisely this internal 
symmetry was also assumed in the original discussion of 
deconfined crticality in 2d quantum magnets with easy 
plane anisotropy^. Generically, either the bosons or the 
vortices are condensed, which implies that either £7(1) 
or time reversal symmetry is broken. However, if the 
critical point separating these states is stable to fluc- 
tuations, then one could tune a single parameter and 
access a deconfined critical point on the surface. It is 
presently unclear if this is true for the theory in Eqn. 
|8l| which is an easy plane NCCPi with a two fold A 
anisotropy term. There is mounting evidence that the 
SU(2) symmetric NCCPi model supports a stable quan- 
tum critical point. Wh ile initial studies wer e div ided 
between continuou d 40 ! 46 ! and weak first ordeJ^EHI re- 
cent studies of quantum models seem to favor continuous 
transitiorPSHSil, Howe ver th e situation is less clear with 
easy plane anisotrop} |4 ° | 54 1 , and the A anisotropy term 
above. The connection to SPT surface states should pro- 
vide additional motivation for further study. 

Thus far we have assumed translation invariance on 
the surface, but in fact only internal symmetries are re- 



quired to define the phase. The presence of surface ran- 
domness that respect internal symmetries will provide 
random variations in the local critical coupling. This 
random energy density term is known to be typically rel- 
evant at a quantum critical pointful, since it requires a 
rather stringent condition to be met, v > 1 for irrele- 
vance in a clean critical pint in d = 2. Here we emphasize 
a crucial difference with the realization of the deconfined 
critical theory in 2d quantum magnets. In that case the 
presence of disorder leads to a random field that cou- 
ples linearly to the VBS order parameter. In the spinon 
representation this is a random monopole insertion term. 
Alternately in the dual vortex representation this is a 
random term that couples to V^+V-^-- This coupling is 
expected to be relevant at the clean deconfined critical 
point and might potentially lead to confinement at the 
resulting disordered fixed point. Thus it is not clear if the 
NCCPi description is a useful one in the presence of dis- 
order. In the present problem however a linear coupling 
to tp2+' l P'2- remains forbidden even in the presence of dis- 
order. The random energy terms though relevant are still 
not expected to lead to confinement by themselves. More 
dangerous potentially are random variations in the cou- 
pling A. The fate of the disordered NCCPi model in the 
presence of this particular kind of randomness remains 
to be investigated. In this context it may be relevant to 
note that even the fate of the 3D fermonic topological 
insulator surface states in the presence of disorder and 
interactions is also not a settled issue. It is currently un- 
clear if one of the symmetries is spontaneously broken in 
the low energy limit. If the symmetries are preserved, a 
critical metal with universal conductance is predicted^. 

As in the other examples a surface state with Z 2 topo- 
logical order that preserves all symmetries is allowed and 
is readily accessed by condensing paired vortices. For 
both Phase 1 and Phase 2 discussed in this subsection 
with symmetry U(l) x Zj, the symmetry properties of 
the corresponding surface topological order is summa- 
rized in Appendix [D] 

Surface Theory of Phase 2 and Vortex Metal: 
A different topological phase is accessed by enlarging the 
symmetry momentarily to [U(l) x U(l)] x Z 2 . Now, we 
can assume that boson of species 1, b\ = e 1 ^ 1 be charged 
under the first t/(l) symmetry and transforms exactly as 
in Eqn|79} Vortices in this boson field transform under 
the remaining U{\) and time reversal symmetry as: 

* 2 -> e"' CTi/2 * 2 : U(l) (82) 
#2 -> UTy9 2 ■ Zl (83) 
Now, the effective theory at the surface is: 

4 = EH- ^)^l 2 + ~/ 2 V + v(|* 2 | 2 )) 

(J 

which has neither monopole insertion nor anisotropy 
terms due to the presence of separate boson species con- 
servation. There is no background flux due to time re- 
versal symmetry. However, breaking this down to the 
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single U(l) symmetry one is allowed the following term 
allowed cos(</>i — fa) since both bosonic fields transform 
the same way under symmetry. Now Vj* = e 1 ^ 1 corre- 
sponds to a monopole insertion operator - thus this is 
a composite operator which in the variables above may 
be written as: ip 2+ fa-V\ + h.c. This breaks down the 
U(l) x U(l) symmetry to a single U(l), and leads to a 
binding of their vortices. Now, if only the time reversal 
symmetry is broken, then a quantized Hall effect results 
of the conserved charge. The discussion closely parallels 
that in Section HvXTl 

However, if U(l) symmetry is broken, the surface is 
a superfluid (properly viewed in the present context as 
an XY ordered state of the spin system), and the vor- 
tices can be shown to be fermonic as in Section HVDI An 
advantage here compared to the bosonic topological insu- 
lator surface is the absence of background boson density 
that implies the fermionic vortices move in zero back- 
ground field, which allows for a sharp definition of their 
statistics in terms of the Berry phase under exchange. 
Moreover, since the vortex density does not break time 
reversal symmetry, generically a finite vortex density will 
be present in the ground state. In the superfluid XY or- 
dered state these vortice will form a vortex solid and their 
statistics is not very important. It is however extremely 
interesting to ask about the result of destroying the XY 
order by melting the vortex solid and proliferating the 
vortices. With Fermi statistics the vortices will form a 
Fermi surface which will be coupled to the non-compact 
U(l) gauge field. The resulting state is a "vortex metal" - 
a compressible metallic phase of bosons with many inter- 
esting properties^. Ref. [28] proposed this vortex metal 
phase as an exotic possibility for a quantum vortex liquid 
in two space dimensions. Here we see that it is a very 
natural outcome for a symmetry preserving phase at the 
surface of this bosonic topological insulator. 

3D Bulk Theory For both Phase 1 and Phase 2, the 
3D topological theories are identical, and only differ in 
the coupling of the conserved charge to the external field: 

£tot = £topo + £em (84) 

C top o = —eBjdcn + <d-^daida 2 (85) 

Under the Z 2 symmetry, an ax, and a 2l — > a 2 i while 
their components change sign. Thus the 'axion' field 
O must be odd under Z 2 so the action as a whole is 
invariant. This allows us to fix O = 0, n, and of course we 
pick the latter value in the topological phase. In general, 
in a 3D topological phase protected by time reversal, both 
fields should transform in the same way. 

Now, Phase 1 has a single charged bosons, fa, 

C-Z asel = ~}-eda 2 dA + ^eBxdA (86) 

47T Z7T 

The bulk theory predicts 9 — bosonic Magnetoelectric 
effect for this phase. 



However, for Phase 2, both bosons are charged so: 

£™ ose2 = ^-e(d ai + da 2 )8A + ^-e(B 1 + B 2 )8A (87) 

and the bulk theory predicts 6 = ir bosonic Magnetoelec- 
tric effect for this phase. 

VII. CANDIDATE PHASES WITH HALF 
QUANTIZED SURFACE THERMAL HALL 
EFFECT 

Thus far we have based our discussion of novel 3D SPT 
phases on the two dimensional K — a x matrix. When 
a conserved charge is present, these phases often lead 
to a quantized magnetoelectric effect, or equivalently, a 
half quantized surface Hall effect. On general grounds 
one may expect additional phases based on the fact that 
thermal transport can also be quantized. In these phases 
chiral modes are expected at the domain walls between 
opposite symmetry breaking regions, that lead to the 
quantized thermal Hall conductance. In this section we 
provide a possible field theoretic description of such a 
phase. Unlike the other examples discussed in this paper 
the issue of whether it can in principle be realized in a 
physical model on a lattice is a bit subtle and is discussed 
at the end of the section. 

Recall that in a 2D system, the combination n xy /T = 

i>t 3h B is quantized, where n xy is the thermal Hall con- 
ductance and T is temperature, in the limit of T — > 0. 
Here vt counts the number of chiral boson modes at the 
edge. For bosons with SRE in d = 2, it is known that the 
quantization takes values vt — 8n that are multiples of 
8 times quantum of thermal conductance. Anything else 
leads to topological order. These states are based on the 
K matrix of the Kitaev Eg stated 
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(88) 



One may utilize this fact to construct the following 
3D SPT phase. Assume time reversal symmetry Z 2 is 
present . Consider the three dimensional theory given by 

C = — V eBjda! + V -if edajdaj (89) 
1=1 I, J 

As long as all the fields aj transform the same way under- 
time reversal symmetry, the coefficient is quantized to 
= 0, 7T. The latter leads to a topological phase. If time 
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reversal symmetry is preserved in the bulk but broken 
on the surface, it is readily seen that each domain has 
thermal Hall conductivity — ±4, and a domain wall 
between opposite domains has the eight chiral edge states 
of the two dimensional theory specified by the K matrix 
in Eqn. [88} 

It is likely that such a state lies beyond cohomology 
classification since we have already identified a phase 
based on K = a x in this symmetry class that exhausts 
the set of states predicted by cohomology theory-^. A 
question that is relevant in this context is whether the 
field theory above can be realized within a lattice model. 
One difference from the other topological phases we have 
described based on K = a x 'FF' term is that in those 
cases a lattice regularization of the field theory can be 
readily envisaged since it involves a Berry phase for the 
product of electric and magnetic fields representing two 
different species of vortices eFiF 2 E^B^ + E^ ■ B\. 
This term is naturally discretized by assuming the corre- 
sponding vector potentials live on the links of the direct 
lattice and the dual lattice. However, the diagonal en- 
tries of the K Es matrix above lead to terms that are not 
obviously compatible with a lattice. Whether this im- 
poses an additional constraint on possible phases is an 
important open question. If indeed the phase described 
above is physically admissible, then it remains to be clar- 
ified if the additional states lead to a Z 2 or a Z exten- 
sion (assuming just time reversal symmetry). We leave 
these questions to future study. In this context it may be 
relevant to note that the analogous free fermion phases 
are topological superconductors in 3D protected by time 
reversal symmetry (Class Dili), which are classified by 
integers. 



VIII. CONCLUSIONS 



describes a putative gapless state from which, on be- 
ing subject to various perturbations, realizes different 
ground states of the surface. It is also relevant to note 
that bosonic analogs of topological superconductors ex- 
ist where domain walls between opposite time reversal 
symmetry breaking regions carry gapless chiral modes. 

There are several open questions for future work. 
Clearly, a central question is whether there are micro- 
scopic models, or perhaps even experimentally relevant 
systems, which could realize these phases. The bulk 
sigma model field theories may provide useful guidance 
in searching for such realizations. One route to accessing 
SPT phases in two dimensions is to start with a frac- 
tionalized phase and confine the fractionalized excita- 
tions. Our analysis suggests that a similar route may also 
be possible in three dimensions by starting with a frac- 
tionalized phase with emergent deconfined U(l) x U(l) 
gauge fields if these are confined by condensational of mu- 
tual dyons (where a monopole of one U(l) gauge field is 
bound to particles that carry gauge charge of the other 
gauge field). Exploring this possibility might also sug- 
gest physical realizations of the 3d SPT phases. A more 
formal question is whether one can push the field theo- 
retic descriptions of this paper to obtain all possible SPT 
phases in 3D, which could shed light on the way in which 
the chiral phases augment the cohomology characteriza- 
tion. The three dimensional 'BF+FF' theories seem a 
convenient tool to capturing bosonic SPT phases. How- 
ever, general constraints on the form of such theories is 
presently unclear. For example, even for bosonic topo- 
logical insulators, does the requirement of lattice regular- 
ization further constrain allowed K matrices that appear 
in the 'FF' term? It would also be interesting to clarify 
the constraints on K matrices when topological order or 
fermions are present in the bulk. 



In summary, we would like to highlight the remarkable 
similarities between free fermion topological insulators, 
and the bosonic interacting topological phases described 
here. In the former case, the surface is gapped only on 
breaking one of the defining symmetries of charge con- 
servation or time reversal symmetry. Then, the resulting 
ordered phase also possesses unusual properties, for ex- 
ample when charge conservation is destroyed by a super- 
conducting surface, the vortices carry a Majorana zero 
mode. Similarly, for the bosonic topological insulator 
with the same symmetries, breaking charge conservation 
at the surface leads to fermonic vortices (albeit with- 
out an attached Majorana zero mode). On the other 
hand, breaking just time reversal symmetry leads for the 
fermionic case to a quantized magneto electric effect of 
8 = t:, whereas for bosonic TIs in the same situation, 
the same response is quantized, but at 9 — 2tt. The 
fully symmetric surface of the fermonic TI, from which 
these conclusions can be readily derived, is a Dirac dis- 
persion of free fermions. We propose that the analog for 
bosons is the deconfined quantum critical action, which 
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Appendix A: Derivation of 3D BF Surface States: 

As a warmup let us recall the derivation of the edge 
states of a Chern Simons theory in 2EPJ We specialize 
to the K = a x Chern Simons Theory: 

£cs = ^-ald v a\ (Al) 
in 

Note, gauge invariance at the surface can be ensured by 
working in the gauge eto = 0. This implies da 1 = so 
a\ = di<j>i. This gives the edge Lagrangian: 

C = ^d x ct>xd T cj>2 (A2) 

leading to the usual Kac-Moody commutation relations. 

The edge dynamics originates from other terms. For 
example, we can add a Maxwell term to the original ac- 
tion (df^al — <9„a^) 2 . The only low derivative term that 
will appear at the edge is from d y a x . Substituting the 
edge field a\, = d x cj) 1 and noticing that the derivative 
perpendicular to the edge (i.e. along y), only picks up 
the confining wavefunction of the edge states we are lead 
to: d y a x oc d x ^ . This gives potential terms 

Note, the pair of fields <f>i, <j) 2 which are canonically 
conjugate Eqn. |A2| is like any regular one dimensional 
Luttinger liquid. The special physics of SPT phases 
arises from the fact that the fields can transform under 
the symmetry in ways that a ID system cannot. For ex- 
ample, in the U(l) protected integer quantum Hall phase 
of bosons, the transformation law of the first nontrivial 
phase is: tf>i — > 4>i + e. i.e both fields transform un- 
der the charge rotation. This leads to protected edge 
states. By analogy, it appears that we should find that 
the surface of a 3+1D SPT phase of bosons is a regular 
two dimensional bosonic system, apart from application 
of symmetries. Indeed, we show below that this is the 
surface state of the 3 + 1 D BF theory. 

Let us begin with the following 3D Lagrangian: 

Cbf = —z — Buvdxdcr (A3) 

Z7T 

Here, a bosonic current has been written as = 
e^ Xa d u B\ (y /2'K, and d A a represents the vortex loops. 

To derive surface properties, again, the non-dynamical 
parts of the Lagrangian implements the constraint: 
tijdidj = and etjkdiBjk — 0. One can solve this to 
obtain a, = di<f> and Bij — e^diOtj. We take the gauge 
Boi = ao = 0. The topological part of the edge La- 
grangian (the edge is taken to be perpendicular to z and 
we use the indices a, b — (x, y)): 

C = ^-d a a b d r <i> (A4) 



Now, for the dynamics, once again one introduces 
'Maxwell' terms in the bulk: (da) 2 and (dB) 2 . Again, 
the ones that survive with low derivatives have d z acting 
on them. The three terms that appear are: 



A - Pi[(d x <t>) 2 + (d y cj)) 2 } + p 2 [(d x a y - d y a x ) 2 } (A5) 
. Thus our boundary Lagrangian is (S — J dxdydr C e ): 

C e = —d a a b d T (j) (A6) 
+Pi[{dAf + (dyej)) 2 } + p 2 [(d x a y - d y a x ) 2 ] 

One interpretation of this Lagrangian is that of a pho- 
ton on a 2D surface where the Gauss law constraint has 
been solved i.e. 2 + 1 D, the Gauss law d x E x + d y E y = 
can be solved by writing E a = e a bdb4>/7r. Then, the term 
that leads to canonical quantization of electric fields: 
E a d T a a is now replaced by ^d T <fi(d x a y — d y a x ), the first 
term in C e . The Hamiltonian of the Maxwell theory: 
pi[E 2 +E 2 ]+p 2 {d x a y — d y a x ) 2 is the term written above 
in the Lagrangian. 

Although Refs. 1201211 follow a rather similar deriva- 
tion, they interpret the theory above as the bosonized 
description of a 2+ ID Dirac fermion. This interesting 
speculation does not appear to be compatible with the 
well known fact that the theory described by Eqn |A6| is 
dual to a 2+1D theory of bosons. Explicitly this can 
be seen as follows. Since e ah d a ab/2n is conjugate to the 
phase (/>, we denote it by n^, = e ab d a ab/2n and use this 
field to write the Hamiltonian of the surface theory as: 



H = 47r 2 p 2 n 2 +pi(V0) 2 (A7) 

This is just the theory of a boson in the two spatial di- 
mensions of the surface (as expected, since we began with 
a bosonic theory). 

It is useful to catalog the connection between the dual 
descriptions. The monopole insertion operator in the sur- 
face electrodynamics is e 1 ^ and actually corresponds to 
the insertion of particles. The other excitations, 'gauge 
charges' of the gauge theory, are the ends of vortices of 
the 3D bulk and are point particles on the surface. Here 
they behave like charges in the 2D electrodynamics, since 
vortices of the <f> field are equivalent to violating Gauss 
law for the electric field. Also, since Curia — > CurlV0 
they correspond to the ends of the three dimensional vor- 
tices. The vortex insertion operator is of course a nonlo- 
cal object, which reflects the fact that one cannot insert a 
gauge charged particle in the bulk without changing the 
gauge fields everywhere. The 2D surface is gapped either 
by monopoles e 1 ^ or by vortex condensation (Higgs mech- 
anism). However symmetry may forbid these leading to 
SPT phases. 
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Appendix B: O periodicity in multicomponent BF 
theory 

In this Appendix we prove the 2ir periodicity of for 
the multicomponent BF theory. The Lagrangian is 



C 



1 -aU t iv\ K d u B I Xti + -^K IJ e tn , XK d tl aid x a J K (Bl) 



2tt " "* 



8tt 2 



Summation over repeated component indices /, J, ... is 
implicit. The crucial second term when expressed in 
terms of the electric fields e 1 and the magnetic fields b 1 
takes the form 



e 

8^2 



5>" 



(2e J . b 1 ) + K IJ 2 (e'.P + e J .b 1 



J>i 



(B2) 

Consider the theory on a closed three manifold such as 
a 3-torus of size L x L x L. Through one cycle, say the 
xy cycle, slowly insert 2nnj magnetic flux of species / at 
This leads to a bulk electric field along the 



a rate 



dt 



z-direction: 



1 d(j>i 
~L~dt 



(B3) 



Next slowly turn on 27rmj flux of b\ in the bulk so that 



bl = 



2'Kmi 
L 2 



(B4) 



The quantum amplitude for these processes is given by 
the term in the action and takes the form 

— e i0 IZf ( K " n i m i+J2i^j K I, nim,i) (B5) 

For some particular pair /, J choose nj — 1, mj — 1, and 
all other np — mj> ~ 0. Then the amplitude simply 
becomes e lBK . If all the elements K IJ are integers it 
follows that is periodic under a 2n shift. 



Appendix C: Other Symmetries: U(l)xi Z2 

Here the U(l) can be interpreted as spin rotation sym- 
metry about z axis, while the Z2 is spin rotn. by 180 
degrees about the x axis. 

Surface Theory: Let bosons b\ = e 1 ^ 1 be charged 
under the U(l) symmetry. Then the phase 0i and con- 
jugate number n\ transform as: 

01 0i +e : U(l) 

ni ->• rix : U(l) 

01 -> -<t>i ■ z 2 

m -> -m : z 2 (ci) 

Now we would like to understand how a vortex in a 
superfluid surface state of this boson field will transform. 



The remaining Z2 symmetry will act on the vortices, 
however, it is readily seen this switches vortices to anti- 
vortices. More formally, the vortex fields are coupled 
minimally to the gauge field 0,2 whose flux is the number 
density n\ — (d x a.2 y — d y a2 X ) /2n. Now since the num- 
ber density changes sign under Z 2 , so does the gauge field 
012 —ol-x- This implies that for the minimal coupling to 



remain invariant, we need 'J 2 
transformation is: 



vpS. In fact the desired 



\&5 



z. 



(C2) 



This may be viewed as the single projective represen- 
tation of U(l) x Z2, where the U (1) may be viewed as the 
gauge U(l) that changes sign under Z 2 . It is readily veri- 
fied that the gauge invariant combinations |-02+| 2 — I ^2- 1 2 
and , 02+V , 2- = e l< ^ 2 both transform nontrivially under 
Z2 verifying that if vortices condense they always break 
the symmetry. A second species of bosons is defined by 
b\ = e 1 ^ 2 , which transforms as: 



02 02 
n 2 -> n 2 



7T : Z 2 
Z 2 



(C3) 



and is neutral under the global U(l). This satisfies 
the intuitive requirements of a topological surface state 
and therefore we conclude a U(l)xZ2 symmetry group 
also leads to Z 2 topological phases. Note however, of the 
two symmetries were in direct product, there would be 
no topological phases. 

Let us write down the field theory for the surface in 
terms of vortices of ^2- They are minimally coupled 
to a vector potential a 2 whose flux is the boson density 

V X Ct\ = Til. 



1 

2k 



f 2 

J 2fiv 



v(|f 2 | 2 )) 



Since the field b\ is charged, monopole insertion operators 
are forbidden, but various anisotropy terms, involving 
four vortex fields are allowed. These and other allowed 
perturbations are readily identified given the symmetry 
transformations above. 

A dual description of the same theory is obtained by 
'fractionalizing' the boson field b\ = ip* + ipi_, where 
^1 = V'l-) ma y be viewed eitheil^ as a Schwinger 

boson representation of b\, or as vortices of b 2 . Now, 
these transform under a projective representation of the 
global symmetry U(l) x Z 2 : 



(C4) 
(C5) 



*i e"^- /2 *i : U(l) 
*i -> cr^i : Z 2 



this is compatible with the transformations in Eqn. |C1| 
We see this implies that vortices in 62 carry half unit of 
global charge at the surface. This will help us fix the 
bulk field theory. 



2G 



3D Bulk Theory Given the characterization of the 
surface states above, we can write down a bulk 3D the- 
ory that reproduces these features. We write down the 
following theory based on K — a x where the conserved 
charge is coupled to an external electromagnetic field A, 
and later justify it: 



Field 


zl 










f 


f 



TABLE IV: Zl 



£tot — £topo H~ ^em 

(C6) 

Ctopo = —e(Bi9ai + B 2 da 2 ) + Q-^daida 2 (C7) 

C em = ^e(da 2 )dA+^-e(B 1 )dA 

Under the Z 2 symmetry, B\ — > —B±, a\ — > — a\, but 
B 2 — > B 2 , a 2 —> a 2 . Thus the 'axion' field 9 is odd un- 
der Z 2 so the action as a whole is invariant. This allows 
us to fix 9 = 0, 7r, the latter value yields the topolog- 
ical phase. Also, only one of the boson species carries 
global U(l) charge which implies that there is no topo- 
logical contribution to the magnetoelectric polarizability 
i.e. = 0. 



Here we discussed two phases (Phase 1 and Phase 2): 
for both there are a pair of e and a pair of m particles 
with different charge assignments as given in Tables |V| 
and ED 
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TABLE V: U(l) x Zl: Phase 1 



Appendix D: Symmetry Transformation of Surface 
States with Topological Order 

For convenience we accumulate in this Appendix the 
properties of the surface state with Z 2 topological or- 
der of the SPT phases with various symmetries. As de- 
scribed in the main paper such a surface topologically 
ordered phase provides a particularly simple perspective 
on why a trivial gapped symmetry preserving surface is 
not allowed. The Z 2 topological order has four distinct 
quasiparticlcs which we will denote 1, e, to, /. The trivial 
quasiparticle sector is described by 1 and consists of all 
local operators. We will take e (for 'electric') and to (for 
'magnetic') to be bosons and / to be a fermion. e, to 
and / are all mutual semions. Below we summarize how 
the physical symmetry is realized for each of the three 
non-trivial quasiparticles for the various phases. In what 
follows q demotes the charge under the global U(l) sym- 
metry. In all cases other possible quasiparticles in the 
same sector is obtained by adding trivial quasiparticles. 

Symmetry U(l) x 
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TABLE III: (7(1) xi Zl 

Symmetry 

Here there are a pair of e particles (denoted e a — 
(ef,e±), a pair of to particles (m a — (to^jTO^)) and a 
single / particle. 

Symmetry (7(1) x Zl 
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TABLE VI: 17(1) x Zl : Phase 2 



As we emphasized in the main paper the realization 
of symmetry is such that these Z 2 topologically ordered 
phases cannot arise in strict two dimensional models with 
local action of the symmetry group. This is readily seen 
from the general i^-matrix classification of 2d time re- 
versal invariant gapped abelian phases in Ref. [M) For 
the case of Z 2 topological order we may, as usual, take 



K = 




(Dl) 



In the notation of Ref. |24l we need to specify the T 
matrix (which describes how time reversal acts on the 
Chern-Simons fields), the charge vector r, and the time 
reversal vector \- The latter describes the action of time 
reversal on the quasiparticle creation operators. With 
the 2x2 K-matrix above, T-reversal invariant insulators 
have 



T 




(D2) 



(The ± sign depends on whether the global (7(1) charge 
is even/odd under time reversal). Further \ = ( x j 0) 
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with x = or 1. The case x = 1 describes a situation in 
which one of the two bosonic quasiparticles is a Kramers 
doublet. In particular in the classification of Ref. [Ml 
it is not possible for both entries of x to be 1, i.e for 
both bosonic quasiparticles (e and m) to simultaneously 
be Kramers doublets. Thus states of the the kind de- 
scribed above for the surface of 3d SPT phases with only 
Z2 symmetry (or by extension U(l) x Zj) are not al- 



lowed in strict 2c? in spin/boson lattice models with on- 
site action of the global symmetry. In passing we note 
that this precludes the possibility that strictly 2d spin 
models have gapped Z2 topological phases where both 
non-trivial bosonic quasiparticles carry spin- 1/2 (i.e are 
spinons) while the fcrmionic quasiparticlc carries no spin. 
The case of U(l) xi may also be analyzed similarly and 
shown to not occur in the if-matrix classification. 
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